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I. INTRODUCTION 


A. GENERAL 


This paper is the result of a Headquarters Marine Corps 
request to explore the use of James-Stein and other 
shrinkage type parameter estimation schemes for the purpose 
of generating manpower loss rates. The very large number of 
cells within the  USMC officer force structure leads to tne 
condition that empirical attrition rates are unstable. This 
problem is compounded by the fact that the cell probabili- 
ties are small. Further difficulties are present because 
some or tne inventory cells are empty for structural reasons 
while others are empty by Chance. Tnerefore, the small 
inventory cells draw especial attention. 

In order to deal with this phenomenon, variance stabi- 
lizing and syametrizing transformations rollowed by versions 
of the James-Stein shrinkage technique were applied to 
selected aggregates of the USMC officer data. The applica- 
tion of these methods to the summary data provided by Navy 
Personnel Research and Development Center (NPRDC) has shown 
much fromise. It has been illustrated within this project 
that improvement can be attained from application or these 
methods over other estimating metnods. 

Since considerable further work is needed to screen our 
choices, a further purpose or this thesis is to document a 
number of supporting materials, including a thorouga discus- 
sion cf the following: 

IENEDackuBousd Of attriticn and promotions for USMC 
officers. 


2. Data types and the immediately available data kase. 
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3. Identifying the plans (consumers) that ٦ 
attrition rates. 
4. Formulas and matnematical risk calculations relating 
to the James-Stein and empirical Bayes estimators. 
Cross validation procedures were utilized to illustrate 
the improvement of shrinkage methods over the use of 


empirical rates. 


B. BACKGROUND 


The forecasts made by Manpower planning models are 
affected by three general factors: existing inventory, 
projected losses, and projected gains. In order to projec 
the inventory into various future time periods it is neces- 
sary to use a realistic system of flow rates. Currently 
those rates are generated largely through arbitrary manual 
processes. The goal of this paper is to exanine ways to 
develcp objective and computerizea loss rates. 

There is an inherent confusion in terminology of losses 
and loss rates in that some are leavers from the sarine 
Corps and some are merely transitioning within the force 
structure. Flows of Officers from particular cells are 
characterized by Military Occupational Specialty ΟΠΠ 
length cr service (LCS), grade, and perhaps some other cross 
Classifying characteristics. When the flows go from one 
cell tc anotner within the Marine Corps, then the flow is 
referred to asa strength loss ror tne former ceil. AD 
ofricer moving from one cell te anotner cell does not mean a 
loss from the Marine Corps to tne civilian labor market. 
Changes such as promotion to a higher pay grade or losses to 
a cell due to an cfzicer increasing from one length of 
service to tne next are other forms of transitico within the 


force Structure: 
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Losses from within the Marine Corps to the civilian 


labor market are the only losses treated in this thesis. 


They can be voluntary or involuntary. Voluntary losses 
occur when officers resign, retire, or are released by 
ShoLlce. Involuntary losses occur due to discnarge, death, 
disability, release from active duty, and retirement. 


Retirement can pe voluntary or involuntary. 

Authors who have studied loss rates and how they apply 
to Manpower planning are emphasizing the importance of 
understanding these rates in order to conduct proper 
forecasting. (Ref. 1) (Ref. 2] 

A dirficulty encountered in doing this project was that 
of obtaining sufficient data points with which to conduct 
Cohort ος census manpower data analysis (Ref. 1]. The 
reason for this difficulty was the state of the Marine Corps 
automated data reporting system prior to 1976. During tne 
period 1970 to 1970, tue Marine Corps was instituting its 


first automated personnei reporting system. Many problens 
were encountered during this period. Therefore, data Iron 
this period is generally unreliable. Data at the Defense 


Manpower Data Center is based on input from Headquarters 
Marine Corps and was found to have similar problens. À 
complicating factor during this period was the Vietnam 
conrlict and the attritions generated by the conflict. 
These prcblems prevented an extended data base from which to 
draw frcm ror the purpose of this paper. The Marine Corps 
manpower reporting system was amended in 1976 and by 1977, 


the system was reporting data with over 95 percent accuracy. 
C. THE MARINE CORPS OFFICER ATTRITION AND PROMOTION 
STRUCTURE 


An analysis of a Marine Corps officer as he moves 


through his career reveals key periods in wnich an officer 


is most likely to ieave the service. These periods coincide 
with promotion points in the service. Under normal circum- 
stances, unrestricted 2nd Lieutenants are promoted to 15% 
Lieutenants after Z4 montns in grade. Ist Lieutenants are 
promoted to Captains in the Fifth year of commissioned 
service, Captains to Majors during the tenth year, Majors to 
Lieutenant Colonels during the sixteenth year of service, 
and Lieutenant Colonels to Colcnels during the twenty second 
year of service. An important influencing factor is the 
fact that retirement benefits are not achieved until the 


twentieth year of total service. 


1. Lieutenants 


Lieutenants normally have a four or five year 
initial term of service depending on their source or commis- 
sion and Military Occupational Specialty (MOS), either air 
OC grona. Aviation Lieutenants are some what special. 
Only a smali number of aviation Lieutenants leave the 
Service. This is due to their longer initial term of 
service, which causes them to reach the grade of Captain 
before the level or natural attrition is more iikely to 
GEC UE. attrition of Aviation Lıeutenants is primarily a 


function of accidents, illness, and disciplinary problens. 


2. Captains 


There is a much greater period of rlexibility for 
Captains given that they are regular officers or reserves 
with an extension of service agreement. Captairs are 
promoted at five years of commissioned service and are 
either passed over or promoted to Major by their eleventh 
year cf service. Thus in a period oc over six years, ап 


individual may leave the Service voluntarily at any point. 


16 


3. Majors 


= m 


"He Tank Of Major is -τπρῖες to characterize. 
Pre-DOPMA (Defense Officer Personnel Management Act)  4ajors 
have a guaranteed length of service of twenty years commis- 
sioned service. Since they are promoted during the tenth 
and eleventh years of commissioned service and are either 
passed over or selected for Lieutenant Colonel in their 
sixteenth to twentieth years of service, they have a period 
of six to ten years in which to make a decision about 
leaving. Post-DOPMA Majors ao not have a guaranteed length 
of service of twenty years but there are provisions wnich 
allow the retention of these officers on active duty based 
on their skills and the findings of a special board at 
Headquarters Marine Corps. By the time officers make the 
rank cf Major most have essentially committed themselves for 
a twenty year term of service. Reaching the twenty year 
retirement point assures one a satisfactory retirement 
program. 

At tne sixteenth year of service, a Major is first 
eligible for selection to Lieutenant Colonel. The selection 
process may be repeated for rour years until the officer is 
either promoted or passed over and retired at twenty years 
or service. During this later period it is very unlikely 
that any Major would leave the service unless he had prior 
enlisted service. That would allow him tc achieve a total 


or twenty years of service earlier. 


4. Lieutenant Colonels 


A ee a coe umm mam سے‎ ee O oe 


Lieutenant Colonels are usually promoted at sixteen 
years of commissioned Service and may continue in service 
Until the twenty sixth year or service. Lieutenant Colonels 
incur a three year obligation upon promotion; therefore, at 
any pcint after their third year in grade upon exceeding 


nineteen years of service they may retire voluntarily. 
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р. ATTRITION OF AVIATORS 


The attrition rate of aviation officers was affected by 
the intiation of Aviation Officer Continuation Pay (AOCP) in 
1981. AOCP provides a bonus per year to aviation orficers 
and this in turn obligates continued service. The progran 
was appiiea to all ranks provided the individual met certain 
active duty rlight status requirements. This action by DOD 
has had its desired effect upon the retention of aviation 
officers.: This accounts, in part, for the reduction in 
attriticn rates for Lieutenant Colonels in 1983 and for 


Majors and Captains in 1982 and 1983. 


E. TERAS 


The following terms will be utilized within this 


project: 
1. Aggregation - collection of historical officer data 
over distinguisning characteristics. For exanple, 


the aggregaticn of all grades within a Military 
Occupational Specialty (MOS). 


2. Attrition - any departure from the Marine Corps py an 
officer: 

3. Attrition rate - idealized ratio of attritions ٣ 
inventory. 

4. End-strength - total number of the Marine Corps 


officers at the end of a specified period of tinme. 

5. Limited Duty Officer (LDO) - professionally qualified 
officers specificaily designated for limited duty 
within certain MOS's. 

6. Unrestricted - all otticers not designated as Limited 
Duty Officers (LDO). 

7. Regular - unrestricted officers not specifically 


assigned to a "Reserve" status. 
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8. keserve - unrestricted officers with a specifically 


assigned "Reserve" status. 


NS 


II. STRUCTURE OF THE DATA BASE 


A. SUMMARY DATA FILE 


The data used in this project was obtained from a 
Summary data file from the Commander, Navy Personnel 
Research and Development Center (NPEDC). The Summary data 
file was generated from two source files; the Headguarters 
Master File (HMF) and the Quarterly Statistical Transaction 
File (STATS). The Headguarters Master File (HMF) was used 
to produce nistorical officer inventories as of the begin- 
ning cf the fiscal year. inventories were generated for 
each fiscal year from 1977 through fiscal year 1983, The 
inventories were identiried by distinct characteristics. 
These characteristics were Military Occupational Specialty 
(MCS), grade and length of service. 


The data in hand does not distinguish between unres- 


tricted officers and limited duty officers. For that 
reason, discussion cf the characteristics that separate 
these two types cannot be used advantageously herein. Such 


discussion is included however in anticipation of recelving 
a more discriminating data base in tne future. 

The Cuarterly Statistical Transaction File was used to 
generate historical lcsses. Within tnis file, lf the TYDE 
Transaction Code indicated a loss, then the Effective Date 
of Action field would specify the year and month of the 
loss. Losses were classifiea into eight categories for each 
fiscal year 1977 tnrcugh 1983. The losses were further 
Classified into distinct elements by MOS, grade and length 


of service. 


B. SUMMARY DATA FORMAT 


The summary data file classified the Marine Corps 
officer inventory into 40 Military Occupation Specialties, 
10 grade levels, 31 lengths of service, and 8 loss catego- 
ries. The data format is defined by Tables I ana Il. The 
column containing the letters A through E refer to the 
structural zero problem discussed in the next section. See 
the material surrounding the ccntents of Table V. 

Throughout the remainder of this project, when reference 
is made to a particular Military Occupational Specialty 
(M03), the data code reference from Table II will be used 
instead of the actual MOS. For example, this project will 
refer to the Utilities MOS as number 06 and not 11. Et 
should also be understood that the two digit MOS identifier 
listed in Table II is strictly the Occupational Field iden- 
tifier in the USMC MOS Manual. This results in MOS, for the 
purposes of this paper, being an aggregation of what is 


usually understood to be an MOS cell. 


С. SIRUCTURAL AND SAHPLING ZEROES 


There are two kinds of zero inventory cells that must be 

dealt with: 

1. Structural zeroes: A cell wnose inventory is always 
zero because certain grades and length of service 
combinations should never appear in the Military 
Occupational Specialty (MOS). For exanple, a 
Lieutenant Colcnel with 4 years of service in any MOS 
or an infantry warrant officer in MOS 03 does not 
exist. 

2. Sampiing zeroes: according to the sample data, there 
is no officer in the particular grade and length of 
service combination rcr a given fiscal year and 
Military Occupational Speciaity. This is not а 


permanent condition, but one tnat occurs by chance. 
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TABLE I 
DATA FORMAT 


| | 
| | 
| | 
| COLUMN DATA | 
| 3-4 MOS GROUP (0... 39) | 
| 5-6 GRADE | 
| πο LENGTH OF SERVICE (0... 30 | 
| Sis BEGINNING OF YEAR INVENTORY | 
| 14-53 LOSS COUNTS (ТҮРЕЗ 9 | 
| | 
i 
| CODE GRADE | 
| 
| 0 WARRANT OFFICER all | 
| 1 CHIEF WARRANT OFFICE CWO-2 | 
< CHIEF WARRANT OFFICER (CWOSS | 
| 3 CHIEF WARRANT OFFICER (CWO0-4 | 
4 SECOND СТВОЛЕ МАТЧА 
| 5 FIRST LIEUTENANT | 
| © CAPTAIN | 
7 MAJOR 
| 8 LiEUTENANT COLONEL 
| B COLONEL | 
| | 
| CODE TYPE LOSS | 
| 1 VOLUN TAR Y RESIGNAIION | 
| 2 VOLUNTARY REITREWTEN? 
| 3 INVOLUNTARY - DEATH | 
| ü INVCLUNTARY - DISCHARGE | 
| 5 INVOLUNTARY ھ‎ ۲۸٦٣ { 
| O RELEASE 0 + + +4 4 17٣ 
| 7 DISABILITY RETIREMENT | 
8 INVOLUNTARY RETIREMENT | 
| 


For tne purpose of clarity in this paper, the ncnstruc- 
tural zero region will be referred to as being within the 
Feasikle Region. This region is defined by the lower and 
upper limits of the length of service for eacn grade. Table 
III identifies the minimum and maximum years of commissioned 


Service for the Marine Corps officer grades for the 
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Recall, 


Lamiteds Duty Officer. 


unrestricted as well as the 


the latter type is not separated in the current data. 
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TABLE III 


MINIMUM AND MAXIMUM YEARS-OF-COMMISSIONED 
SERVICE FOR OFFICER GRADES 


Grade Unrestrictoq LDO 
IN YCS MAX YCS MIN το MAX YCS 

WO1 1 4 

CWO2 ü 9 

CRO 3 8 13 

CWHOŅH 12 24 

21Τ Ί y = Е 
LER 1 8 1 8 
CAPT 4 13 2 12 
MAJ 9 20 5 18 
ITCOL 14 23 8 23 
COL 21 30 - - 


ΕΕ. rT 
l ! ! ! 
yt 
ye eS EY eee | 


In order to get the officer's length of Service one must 
add the ofticer's prior enlistment period to the officer's 
years of commissioned service. Since prospective warrant 
officers are reguired to have a minimum of five years of 
service prior to consideration by the officer selection 
board, the logical lower limit for length of service for 
warrant officers is five years. The LOS then would include 
the total of this service and the respective years of 
commissioned service. Additionally, warrant officers are 
required to nave 10 years of service prior to consideration 
to pecome Limited Duty Officers. This affects the lower 
limit of the lengta of service ror Limited Duty Officers. 
Tae logical lower limit for length of service for  LDO's 
would be 10 years. 

The summary data provided for this project did not 
distinguish between unrestricted and Limited Duty Officers. 
Additionally, the length of service was the characteristic 
provided instead of the years of commissioned service for 


each of the officer grades. Table IV identiries the lower 
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and upper limits of the Feasible Region of the length of 


service for Marine Corps officer grades utilized in this 


project. 
BM S O. u or 
| | 
| TABLE IV | 
| MINIMUM AND MAXIMUM LENGTH OF SERVICE | 
| FOR OFFICER GRADES 
| Grade MIN LOS MAX LOS | 
WO 1 9 20 
CW O2 8 30 | 
| CWO3 12 30 | 
| CW O4 16 30 | 
2t 1 20 | 
TET 1 20 | 
CAPT E 30 
MAJ J 30 
CI OL 14 30 
COL 21 30 | 
je | 
———-——— ————————-—————-—-———-_—-_—-———.Ї 


Note the maximum limit іп the majority oí grades can 
feasitly be 30 years. This was substantiated by the actual 
data wherein the maximum length of service for grades CWO3 
through Colonel was 30 years. The maximum length of service 
for C1's, 2Li's and 12,745 was limited to 20 years because 
ot the very small number of actual officers with 21 to 30 
years of service within these grades. 

One method of ccnsidering the size or the feasible 
region 15 to create a multidimensional array using all 
combinations of different characteristics within the data. 
If all characteristics were considered tor all seven fiscal 
years, the size or the multidimensional matrix wculd be 
1x40x31x10x3. This would represent 7 years, 40 MOS's, 31 
bes 's, 10 grades and 8 loss types. This would create 
' 694,400 different cells that could be dedicated to loss 
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2600 7۷ However, to bring the dimensicn to a more manage- 
able level, consider one year ort data for each loss type 
having a cell count of 40x31x10. Therefore, for each loss 
type within a fiscal year there 15 а total of 12,400 cells. 
Many cf these cells are infeasible and are designated those 
witn structural zeroes. 

The region ot structural zeroes 1s the area outside of 
the minimum and maximum length of service for the officer 
grades as listed in Table IV. But it 15 further articulated 
because of additional constraints imposed by the various MOS 
classifications. This is largely because of the Suitability 
of certain MOS's being assigned to unrestricted, warrant and 
Limited Duty Officers. The particular MOS*s аге ideéentiried 
for illustration purpcses by letter in Taole II according to 
the general categorization listed in Table V. This general 
MOS categorization specifies by category the inclusive 
grades within certain “OS's, the number of MOS's with the 
specified grade structure, the number of structural zeroes 
per OS, and the total structural Zeroes within each 
category. The total number of structural zeroes is 6149. 
Therefore, the number or cells within the feasible region is 


6251 for each loss type for a given fiscal year. 


D. CENTRAL DATA 


The tollowing iS an observation made by the author 
regarding the summary data utilized in tnis project. This 
data, which was provided by Navy Personnel Research and 
Development Center (NPRDC), was classified as central data 
according to Bartholcuew [ Ref. 1: DS Bartholomew 
derines central wastage (attrition) rates as follows: 

1. The central wastage (attrition) rate is the number of 
leavers during the period who were in this class when 
they left divided by the average number in this class 


during the period. 
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| | 
i TABLE V | 
| STRUCTURAL ZERCES CATEGORIES | 
| SEE: Total | 
| Grades Number Zeroes zeroes | 
| Category Within MOS ӨТ MOS per MOS per s 
| A WOT. s- CICOL 23 129 2367 | 
| В noie C ROS LDO 8 159 1272 | 
С ΠΝ LICOL 3 202 606 | 

| D 521222564 5 231 1185 | 
| Β WOT. COL 1 119 119 i 
| TCTAL 40 6149 | 
um { 
E СЕЕ 3 


A problem arose on several occasions when the data was 
disaggregated to a level where the inventory was very small. 
For example, when examining the inventory in a particular 
fiscal year, the inventory was zero tor a Length of Service 
(LOS) and Military Occupational Specialty (MOS) combination. 
Examining the inventory in the next fiscal year for the sane 
LOS and MOS combination also was zero. However, the problen 
arises when the numper of leavers iS egual to or greater 
than cne. The average of the inventories for the two fiscal 
years is zero. USing this result in the estimation of the 
attrition rate would бе ambiguous. The zero in the denomi- 
nator of the ratio of leavers to average inventory results 
in an undefined expression. 

For tne purpose of removing tnis ambiguity in the data, 
the following policy was adopted to define the central 
inventory number for the officer force at disaggregated 
levels: For any cell cr collection of cells 

1. Iet, t = 1...6, refer to tre years 1977...1982. 
2. Let, Y(t) = Number of l1csses in year t. 


3. Let, INV(t) = Inventory in the beginning of year t. 
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4. Let, N(t) = Maximum of Y(t) and the average inventory 
in year t, using the beginning inventory in year t 
and t+1 and ccmputing their average: 
(INV(t) +INV(t+1)) /2. 

5. Let, m= Sum cf Y(t) divided by sum of N(t), (both 
suns over t) represents the empirical central 
attrition rates for the particular cell or aggregate. 

The data bank that will ultimately be utilized in the 
parent project at NPRIC, of which this project is a subset, 
will contain the appropriate inventory and attrition data so 
that tne above modifications will not be required. It is 
currently planned, ty NPRDC, to produce the historical 
officer inventory as of the beginning of a month and tne 
losses will be identified by the month in which the officer 
leaves tbe Marine Corps. This will provide the elements for 
a more accurate estimation of tne attrition rate on the 


disaggregate level. 
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A. OVERVIEW 


In order to gain an appreciation for the Complexity of 
the Manpower Planning process, the author during his experi- 
ence tour at Headquarters Marine Corps had the opportunity 
to contact key individuals responsible for the execution of 
tne mcdels currently utilized within the Manpower depart- 
ment. The methodology used to identify the Current system 
included interviews with the action officers of the Officers 
Plans Section (Code MPP-30), briefings and interviews with 
representatives from Allocations Section of Manpower Control 
Branch (Code MPC-20), Officer Assignment Branch (Code MMOA), 
Manpower Plans, Programs and Budget Section (Code MPP-40), 
Manpower Management  Informaticn Systems Branch (Code MPI), 
Assistant Deputy Chief of Staff for Manpower (Code M), and 
examination of existing manuals, Requirements statement for 
Orficer Planning Systen, and Statement of Work for the 
Develcpment of the Officer Inventory and Manpower Flow Data 
Base Systen. 

The goal of the officer planning process is to shape the 
officer force structure over the 7 planning years (current, 
budget, and five Program Objective Memorandum (POM) years). 
This fcıce structure is affected by three influencing 
factors: existing inventory, projected losses, and expected 
gains. The Officer Plans Section (Code MPP-30) supports the 
Manpower Department at Headquarters Marine Corps by 
preparing plans that address each ot these factors. 

The plans which input the target inventory calculations 
are the Promotion Plan, the Lateral Move and Directed 
Lateral Move Plans, the Inventory Projection, and the 


Manpower Plan. 
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Losses in tne otficer force are tàe result of retire- 
ments, resignations, releases, discharges and other factors 
(e-g., death). Among these factors, releases and aandatory 
(statutory) retirements addressed in the Oficer Inventory 
Control Plan, can te estimatea with a high degree of 
certainty. The remaining loss factors are estimations based 
on fcrecasted attriticn rates utilizing historical data. 
All of the factors contributing to losses in the officer 
inventory are addressed in the Inventory Projection Loss 
Plan. 

Gains in officer force are tne result of accessions and 
augmentation. Àccessions are identified in the Accession 
plan ana augmentations are “¡identified in the Augmentation 
Plan. The number or field grade officers are directly 
influenced by statutory limitations specified in Title 10, 
ο ο. «Codes These limitations affect both of these plans. 
In the planning process, in order to plan for an adequate 
number of field grade officers in the future, one must 
currently access a sufficient number of second lieutenants. 

The force structure is also affected by the mixture of 
skills of the officer inventory. The plans used to support 
the required nix and levels of skills are the T3S/40S 
Distribution Plan, tne Entry Level Training Plan, tne 
Laterai Mover Plan, the Continuation Plan, and the Promotion 
Plan. 


В. SCURCES OF INFORMATION 


Input from the following external sources provide infor- 
nation required by the Officer Plans Section to create the 
plans identified above. 

1. Statutory limitations specified In Тісіе πο g. 52 
Code. 
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2. Desired inventcry force structure as enumerated іп 
the Grade Adjusted Recapitulatior (GAR). 

3. Feedback data from agencies which exəcute the plan 
(for example, actual monthly promotions). 

4. Historical data from JUMPS/MMS on officer service 
characteristics and any changes in the officer 
invertory over time. 

5. Current data оп the existing officer inventory 
utilizing the Cfficer Slate file andthe Lineal List 
file. 


С. FLOW OF INFORMATION 


The planning process consists of three major elements 

comprised of several frocesses. These elements are: 
1. Preparation for planning. 
2. Plan development. 
3. Force structure analysis. 

The preparation process includes development of the 
targeted force structure to include distribution of officers 
by skill, sex, unrestricted/restricted categories, and 
grades based on historical fiow rates and USHC policy objec- 
tives. The requirements definition process establishes the 
authorized number of officers in eacn Military Occupational 
Specialty (MOS) and grade for the rive planning years and 
determines tne target force. The programmed reqguirement is 
Specified in the Numerically Adjusted Recapitulation (NAR). 

The planning process also includes the determination of 
orficer attrition rates. Currently the attrition rate 


determination method is accomplished in conjunction with the 


Officer Promotion Planning Process. This model divides the 
attrition rates into two categories: Statutory and non- 
statutory. The statutory attrition Trate results from USMC 


and Гор policies and is applied directly to the officer 
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inventory. The non-statutory rate is derived from a tine 
series analysis of officer inventory. This rate is baseå on 
historical data and assumes ccntinuation of present trends. 
Currently there is no attempt to provide predictive rates 
based on econometric techniques. However, such an approach 
is being researched by the Navy Personnel Research and 
Development Center ror use py the Manpower Department. 

The development process includes projection of the 
officer inventory and generation of various plans. The yoal 
of the inventory projection process is to provide a break- 
down by month of the number of orficers in each grade, 
Military Occupational Specialty (MOS), and other officer 
attricutes (for example, restricted, years of commissioned 
service, sex, etc.). This breakdown is for each ot the 
planning years. Currently the inventory projection is 
detertrined in one of several ways, depending on the plan 
being developed. The Promotion Planning Process projection 
is utilized if the level of aggregation is adeguate ror the 
plan. Otherwise, manual calculations are used to aggregate 
trom tne promotion planning process inventory projection or 
to generate a new projection. 

The plan generaticn process 1S complicated by the inter- 
dependence of tne plans and tae timing of their production. 
Table VI identizies the plan generation process including 
each plan, the frequency of the generation and update of 
each plan, the users, and the sources of data. 

The force structure analysis process involves уег1ї1са- 
tion cf the generated plans and evaluation or their quality. 
Comparison is made cf the results of the planning process 
with the USMC objective force specified by the Grade 
Adjusted Recapitulation (GAR) and the promotior flcw points 
specified ру Defense Officer Personnel Management Act 
(DOPMA). 
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| TABLE VI 
| OFFICER PLAN GENERATION SUMMARY 
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| 
| 

SOURCES | 
PLAN GENERATED UPDATED USED BY OF DATA | 
| 
Prcmotion plan Annual 7ت5‎ ۲ ٦77۰۰٠٠٢۶ ۲۹ 3۰ 5 7۰39۰1 20301 | 
MPP-30 б | 
RES | 
MMPR I 
Augmentation Annual Monthly  MPP-30 1522 35 5; | 
plan MMOA 11 | 
ORB | 
Inventory Annual Annual ЯРР-30 12/39/74: | 
Prcjection а | 
Inventory Proj. Annuaı Annual MPP-30 152.1 | 
Loss Plan | 
Continuation Annual Annual MPP-30 l5 || 
Plan | 
Accession Plan Annual Monthly MPP-30 16 2503,4, | 
MMRD 6,8 | 
Entry Level Annual Monthly MMOA 1,2,4 | 
Training Plan 185 | 
1 | 
i 
۲۳۰7/1903 |. Annuai Mortnly  MMOA ац | 
Distribution 155 
Plan T | 
Lateral Move Annual Monthly  MPP-30 225550731 
plan MMOA | 
Manpower Plan Annual Monthly MPP-30 lilerari; 
МРР-40 5٤ء‎ Θ..ῃ | 
ο icer control Annual Monthiy  HPP-30 12224975 | 
an 
| 
{ 
| 
| 
Source of Data Codes | 
1. HMF | 7. JUMPS/MMS | 
2. Attrition Rates 8. WO Accessions | 
3. ФОРМА 9. Retirement/Separation | 
4. GAR 10. Lineal List | 
9. Budget | ٦٦٦02217660۲7518 
6. Candidate list | 
. 0 ση 
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D. PLANS 


Promotion plan is generated by tne use of the Promotion 
Planning Process model. The model's output is a statistical 
projection of the officer inventory over a period of up to 
10 years. It also provides several built in promotion 
strategles. 

Augaentation Plan is generated to provide the Officer 
Retention Board (ORB) with information on the vacancies 
available in each fiscal year of the planning process 
regarding officers in a category (i.e., Aviation, Ground, 
and Naval Flight Ofricer) with a certain number of years of 
commissicned service. 

Inventory Projection Plan is generated for a five year 
projection and is specified ror each year by grade, years of 
commissioned service, and category (unrestricted or Limited 
Duty Cf ficer- (LDO)) 

Inventory Projection Loss Elan contaıns Jata on expected 
losses from the inventory projection process. The plan is 
utilized to ındıcate which year of commissioned service the 
losses are most likely to occur. 

Officer Inventory Control Plan provides information on 
the number of reserve officers whose contracts expires 
during the riscal year. τι 15 influenced by the 
Augmentation Plan and provides guidance on the number of the 
orficers who may be retained. 

Ccntinuation Plan specifies the number of officers tnat 
пау Ге granted an extension of active duty depending upon 
their Military Occupational Specialty (MOS). 

Accession Plan identifies the desired number of incoming 
potential officers and receives as input the number of 
vacancies for each year of the planning process. since 
accessions ınto the USMC are ty officer candidates wno have 


usually Signed contracts several years prior to graduation, 
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the pcoi of accessions for the first three years are gener- 
ally well known. The goals for the remaining planning years 
are established on prcjected officer losses and or rudgeted 
end strengtnhs. 

Entry Level Training Plan and the TBS/MOS Distribution 
Plan ensure that the entry level officers have the mix of 
Military Occupational Specialties that is needed to fully 
support the torce structure. 

Lateral Move Plan provides information on the shortage 
or overage in each Military Occupational Specialty by 
year-group. 

Manpower Plan is a summary of gains and losses in each 


grade. 


IV. MARINE CORPS OFFICER ATTRITION MODEL 


—— eee eC ee MEER er < 


А. DEFINITION 


A definition for a manpower model offered by Bartholomew 


lS: 


A manpower model is a mathematical description of how 
Change taxes place in the system. First of all th= 
requires the specification of any constraints under 
waich the system operates . . . . 3ے‎ ٦ a model 
must SS the mechanism which generates flows. Some 
flows such as prcmotion or demotion are under the 
direct cortrol of Danagement . . . . Other flows, suca 
as Ὃ ο wastage, are not under direct management 
control and assumptions abcut their future levels is 
likely to be based on a blend of historical datan 
management judgement. (Ref. 1: p. 


B. ASSUMPTIONS 


Assumptions regarding flows ot personnel within a system 
can be classified in various ways. If a manpower policy 
dictates that a specific percentage ог individuals in a 
particular category would leave in a given time period then 
the assumption about the rlow would be deterministic. Given 
a level or inventory in the category under the deterministic 
assumpticn, then there is no uncertainty about how many 
individuals will actually leave. A deterministic model 
contains no random variables and there is a unique set of 
model output data for a given set of inputs. ([ Ref. ὃς р. 3] 

However, if each individual in a category had a certain 
probability of leaving in a given time period, then the 
assumption about the flow would be stochastic. Given a 
level of inventory in tne category under the stochastic 


assumption, then a prediction of the number of individuals 
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leaving woula not be precise. A stochastic model contains 
one or more random variables. The output of a stochastic 
model is random and thus only estimates the true character- 
istics of the system [Ref. 3: p. 3]. The output of a 
stochastic model is a probability distribution. If individ- 
uals within a system behave independently regarding their 
decision of leaving, the actual number of leavers has a 
binomial distribution whose average ıs a certain percentage 
of the inventory within a cateyory [Ref. 1: р» Т 
Voluntary attrition from the Marine Corps is a stochastic 
flow since it results from a number of more or less indepen- 


dent individual decisions. 


C. ATTRITION RATES 


A fundamental role in manpower analysis is played by 
attrition rates. The time period can be any duration. The 
most ccmmonly used periods are montas and years. The use of 
tne central attrıtion rate when the Categories are defined 
witn respect to length or service can be considered as esti- 
Dating the constant rate of leaving or separation during a 
year. Using this method, the standard error of the central 
attriticn rate can be estimated ir some assumptions are made 
about the attrition process: 

1. The number in the category throughout the year was 
constant and equal to the average number in tnis 
class during the year. 


2. Let, na = the central attrition rate. 


3. Let, 5% = Average numter ina category during the 
year. 
ü. Let, L* = Numter of leavers during the year who were 


in this category wnen they lert. 
ο Then, m = BDF / 98. 
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6. Each individual in the category is subject to the 
probability  (mdx) of leaving in each interval 
(X,x+dx). LR OEP 

Under these assumptions the number of losses would be 
distributed Poisson (S*m). The estimated standard deviation 


or the central attrition rate would be: [Ref. 1: p. 25] 


sd (m) = m / sqrt (L*) 


D. MARINE CORPS OFFICER CENTRAL ATTRITION RATES 


Scme empirical attrition rates are provided for the 
pur pose of familiarization. First are macro (large aggre- 
gate) examples and proceed progressively to more refined 
categcrized examples. An example of the estimated central 
attrition rates and their estimated standard deviatıons are 
provided in Table VII. The average inventory over the seven 
year period is alsc provided in Table VII in order to 
illustrate the number of officers within each MOS. 

The rates in Table VII are an aggreyatıon of all grades 
within a Military Occupational Specialty (MOS) and all loss 
types. The time period includes all seven years cf inven- 
tory counts and the losses for friscal years 1977 to 1982. 


The central attrition rate was calculated as follows: 


1. Let, t = 1...6, covering years 19/75 SB 

2. Let, Y(t) = Number of losses in year t. 

3. Let, INV(t) = Inventory in the beginning of year +. 
4. Let, N(t) = Average inventory in year t, using the 


beginning inventory in year t and t+1 and computing 
their average: (INV (t) +INV (t+1))/2. 
5. Let, m = Sum of all Y(t) divided by sum of all N (DRE 
The central attrition for each Military Occupational 
Specialty is iliustrated for the iength of service period 
from Q to 30 in Appendix А, in Figures A.) Cosme These 


figures illustrate tne aggregate central attrition rate ror 
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TABLE VII 


BATES AND STANDARD DEVIATIONS 


ATTRITION 


CENTRAL 


RATES 


MOS 
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„001 


„ 102 


18780 


TOTAL 
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1977 το 


the six years from 
These figures demonstrate the breakpoints and discon- 


types of 195565 over 


all eight 


92. 


tinuities within each MOS over the distribution of length of 


One can very easily note the increased rates wnicn 


service. 
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occur at certain years of length of service. The increased 
rates frequently occur at the four, 20, and 30 years of 
length of service. Ihe central attrition rate for each 4OS 
is illustrated for the grades from Warrant Officer to 
Colonel (labelled 0...9 on the abscissa) in Appendix A, іп 
Figures A.11 to A.20. These figures also illustrate the 
aggregate central attrition rate for all eight types of 
losses over the same six years. 

The central attrition rate for each fiscal year aggre- 
gated cver all categories is provided in Table VIII and 
illustrated in Figure 4.1. The trend is stable for the 
rirst three years and a decreasing trend is evidenced in the 


last three years. 


(UR Puen RS ee ee Ee —— D ———— атады 2) 


TABLE VIII 
USMC OFFICER LOSS BATES PER FISCAL YEAR 


1977 1978 1978 1980 19837 1982 
„1133 „1137 21136 10727 сон 


x 
| 
| 
LS 


| 9 


Table IX is provided in order to gain an appreciation of 
the number of losses within each Loss type. This table 
includes the number of losses for each fiscal year 1977 
through 1983 aggregated over all grades, LOS, and MOS. This 
table also identifies the loss types. One can note the 


large number of losses throughout the seven years in types 


two, four aid Six: Additionally, after reaching a total 
higa τα πο, the general trend is a decreasing number of 
losses. 
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USMC OFFICER LOSS RATES PER FISCAL YEAR 
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Figure 4.1 USMC Officer Loss Rates Per Fiscal Year. 


The central attrition rates tor this level or aggregation 
tor fiscal years 1977 through 1982 are iliustrated in Table 
X. 

The certral attritrior, rates For each Military 
Occupational Specialty (HOS) aná ror each Liscal year 1977 
۲--031 1982 іс ргочійес іп Table ХІ. These rates are 
aggregated over all lcss types. Tnese rates are graphically 
illustrated in Appendix A, in Figures A.21 to A.28. One can 


observe tne chanye of rates over time ror eacn MOS. ШЕ сап 
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TABLE IX 
LOSSES PER PISCAL YEAR BY LOSS TYPE 


LOSS | 
| TYPE 1977 1978 1979 1980 1981 1982 1995 ΟΙ | 
| 

D 118 90 09 40 24 37 35 453 | 
| T2 675 701 681 594 531 528 622 ۷۶۷۵ 

ЕЗ 32 42 38 39 34 22 27 234 
| ІЧ 364 515 536 572 532 8ЯЗІ 317) 7 | 
| 15 Ί 1 0 0 1 0 0 3 | 

8 80; 772 750 595 522 406 436 8375 
| 17 15 27 17 2 3 2 10 76 | 
| 18 2 50 327 5 43 30 4 1 290 | 
| TOTAL 2151 2158 2170 1877 1690 1352 MLS TZ | 

| 
| 
| | 
| LOSS TYPES | 
| Lis Voluntary Resignation | 
| L2: Voluntary Retirement 
| L3: Involuntary = Death | 
i 14: Involuntary Discharge | 
| L5: involuntary Pisapri oy 
| Lo: Release trom Active Duty | 
L7: Disability Retirement 

| 8: Involuntary ketirement і 
ie | 


be noted with interest tnat the MOS's with larger invento- 
ries (refer to Table VII for the average inventory witnin 
each MOS) do not change as much Over time as the MOS's with 
smaller inventories. For example, the infantry MOS 03 has 
the secona largest inventory ccunt. The range of the rates 
over time fluctuate from a low of 0.070 to a high of 0.106 
giving a range of 0.036. The logistic MOS 04, on the other 
nand, has a much smaluer inventory count and the rates over 
time fluctuate from a low of 0.049 to a high of 0.22, giving 
a range cf 0.171. In cther words, the total inventory count 
in a category has a definite efrect on the change of rates 


Over. time, 


TABLE X 
CENTRAL ATTRITION RATE PER FISCAL YEAR BY LOSS TYPE 


] 


LOSS 

TIPE 1977 1978 1979 1980 1931 1982 
L1 .006 .005 .006 .002 .001 .002 
L2 .030 5.036 2035 ив .028 
L3 ΤΟ m. 002 002 72002772901 
L4 „ОТО 027” „028 050 5 323 
15 „000 .000 .000 .000 .000 .000 
L6 .,0u7 .0450 .039 .031 .027 .022 
L7 200058500125 001 090 20290 -000 
L8 2005800350029 2902. 720027 002 


AE E |‏ سے ana- Bre‏ سے SIE‏ ےھ سے 


The change in the rates are also affected Гу the level 
of disaggregation. The centrali attrition rates tor the 
major Military Occupational Specialty groups (Avıation, 
Compat Support, Ground Combat and total) are listed in Table 
ХІТ. These major groups are illustratea in Appendix A, in 
٣۰۶2625 A.297 A.30, A.31, and A.32. The rates denmcnstrate 
the aggregation over all types of losses. For brevity the 
illustrations were linited to the grades of Second 
Lieutenant to Lieutenant Colonel. 

A further disaggregation is listed in Table XIII and 
illustrated in Figure 4.2. These rates demcnstrate the 
aggregation within the Combat Support group over all lengtns 


ot service (LOS) for the voluntary retirements (type 2 loss) 


for the grades of First Lieutenant, Captain, Major and 
Lleutenant Colonel. A common ooservation within this loss 
type is the increased rate as tne grade increases. This 


phencsenon also corresponds to the increased rate or 
attrition in the larger LOS categories. 

Another level of disaggregation of the Ccmbat Support 
group is listed in Table XIV and iilustrated in Figure 4.3. 


These rates demonstrate the aggregation within tne Motor 
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TABLE XI 
RATES FOR FISCAL YEARS BY MOS 


ATTRITION 


CENTRAL 


1977 1578 1979 1980 1981 1982 


MCS 
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voluntary retire- 


for tne 


LCS 


for the 
and Lieutenant Colonel. 


cver all 


Transportation MOS 


First Lieutenant, 


grades of 


2 loss) 


nents (type 


Captain, 


Ihis illustrates 


Major; 


voluntary att m 


distribution resulting from 


the irregular 


even with aggregation cver LOS, 


τση 


The number of losses, 


The actual losses іп tn1s table range 


15 relativeiy small. 


from zero to seven. 
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A rurther disaggregation tc include each LOS 


trated), 


ہے 





TABLE XIII 


COMBAT SUPPORT GROUP VOL. RET. NIS TYPE 2) 
CENTRAL ATTRITION RATES 


0تت 


GRADE 1977 1978 1979 1980 1981 1982 | 

T -004 .005 .000 Ὁ 4002 7000 | 

CAPT „031 .017 .023 .017 7027 ЖС | 

MAJ 4103 .105 .092 сСТІЗ ООО ۷ 

LTCOL .122 .178 „11% Ια ΤΡ, ἱ 
elements witnin a cell. This iS particularly evident when 


disaggregation 1S within a Sparsely populated MOS and the 
desired rate is at a particular grade, LOS, and loss type 


combination. 
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Central Attrition Rates. 
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Figure 4.3 Motor Transport Vol. Ret. (Loss Type 2) 


Central Attrition Rates. 
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V. JAMES-STEIN TECHNIQUES AND VALIDATION PROCEDURES 


A. BACKGROUND 


The first tasic process in statistics is the simple act 
of counting and the second most basic process is averaging. 
One use ot observed averages is to predict or estimate unob- 
served quantities. For example, a baseball player who gets 
eignt hits in 25 times at bat is said to have a «320 batting 
average, wnich is the ratio of 8 to 25. In computing 071 
statistic ar estimate is formed of the player's true batting 
ability in terms of his observed average rate of success at 
the plate. A common reply to the inquiry regarding tnis 
baseball player's next 100 times to bat would be that he 
would probably get 32 more hits. In traditional statistical 
theory, no other estimation rule is uniformly better than 
the otservei average. The attrition in cells is modeled as 
a set of Bernoulli trials whose number is the inventory of 
that cell. The process of choosing leavers is viewed as 
independent from individual to individual. A convient 
anology might be that of a baseball players batting average 
in which each time at bat is viewed as an independent 
Bernculli trial. 

Using tne above Bernoulli (Binomial) model, the rates 
(or batting averages) are the Maximum Likelihood estimators 
of the Bernoulli parameter. The maximum lixelihood proce- 
dure examines the likelihood function of the sample values 
and takes as the estimates of tne unknown parameters those 


values that maximize the likelihood runction [{Ref. 4: р. 


These estimates are cptimal (in particular admis-‏ ٭( 63د 
sible) if we are talking atout one cell (ог one ball‏ 
player). But it is no longer true if we are dealing witn‏ 


three or more celis (cr ball players). 
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Charles Stein (Ref. 5] showed in 1955 that it is 
possible to make a uniform improvement on the maximum like- 
lihood estimator (MLE) in terms of total squared error risk 
when estimating several means from independent normal popu- 
lations. In 1901, Janes and Stein [Ref. 6] presented to the 
Berkeley Symposium an estimator which shrinks the maximum 
likelihocd estimator toward the origin, and can improve 
(1.86. ¿Lower risk) on the maximum likelihood estimator quite 
substantially provided there are three ог more cell. To 
continue the baseball example, if one considers three or 
more Laseball players and desires to predict future batting 
averages for each of them, then there is a procedure that is 
better than simply extrapolating from the separate averages. 
Efron and Morris (Ref. 7] demonstrated an application of the 
Stein rule and its generalizations applied to predict base- 
ball averages, to estimate toxomosis prevalence rates, and 
to estimate the exact size ot Fearson's chi-square test with 
results from a computer simulation. Tae result of these 
applications was a mean square error of the Stein rule being 
less than naalf or that when one uses the individual 
empirical averages. 

The first step in applying Stein's method is tọ deter- 
mine the average of the averages. This grand average or 
grand mean must lie between 0 and 1. The essential process 
in Stein's method is the "shrinking" of all the individual 
averages toward this grand mean. When the Stein's method is 
applied to attrition rates, the results are as follows. 

1. If an aggregated rate is nigher than the grand mean 
then it will be reduced. 
2. If ar aggregated rate is lower than the grand mean 
then it will be increased. 
Figure 5.1 is a simple illustration of this "shrinking" 
concept utilizing fictitious values ror the empirical 
values. In this simple illustration the erfective shrinkage 


1s 25 percent. 
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B. JAMES-STEIN ESTIMATOR 


The James-Stein estimator for each desired cell is found 
through the following equation. 
1. Let Y* = the grand mean. 
2. Let m = sample attrition rate tor each cell. 
3. Then, p = Y* + c(m-Y*), for a specialized value of 
"ch, is the James-Stein estimator. 

See Appendix В, for the James-Stein estimation algo- 
rithm. See Appendix C, for more mathematical details on the 
equation used for the James-Stein estimator. The amount 
(n-Y*) is the amount by which tae sample attrition rate 
differs from tne grand mean of the desired cells. The equa- 
tion thus states that the James-Stein estimator "p" differs 


from tne grand mean by tnis same quantity (m-Y*) multiplied 


by a coefficient, o Ihe coefficient "c" is the 
"Shrinking factor" used in the James-Stein estimation 
process. If "c" was equal to 1, tnen the resultant 


James-Stein estimator for a given cell would be identical to 
the celi's sample attrition rate. In otber words, р-п. 
Stein's theorem (Ref. 6] states tnat the shrinking factor is 
always less than 1. Its actual value is determined by tne 
collection of all observed attrition rates. 
As an example, if the shrinking factor "c" was equal to 
„3, then each attriticn rate will shrink about 70 percent of 
the distance to the grand mean. Three examples follow: 
1. Assure the sample attrition rate was .10, the grand 
mean was u, апа: πεν equal to 255% Тһеп 
ο. οσοι οἱ. 102.05). Ihe result is «065. Ey Stein's 
theorem, the attrition rate 1S best estimated by .065 
rather than the historical sample attrition rate. 
Shrinkage was dcwn toward the grand mean. 
2. Assume the sample attrition rate was .50, the grand 


mean was ο. and "ey equal to en Then 


Οἱ. ΠΣ The result rs NEM The esti- 
nator is smaller than the sample attrition rate. 
Shrinkage was dcwn toward the grand mean. 

3. Assume tne sample attrition rate was .05, the grand 
mean was . 10, and wee" equal to v Then 
p=z10+. 7 (2?05= 7۷ The result 77572065: Shrinkage 
was up toward the grand mean. 

There are several expressions for the James-Stein esti- 
nator, Lut all of them have in common a shrinking factor 
Mom. without exhibiting any of the various formulas for 
"с, the following ıs a description of their general 
behavior. See Table XV for tne James-Stein formula utilized 
in this Project, See Appendix C for further background and 
formulas. 

Let "К" be the number of unknown means; let sigma be the 
within population standard deviation; and let SSB be the sun 
of square о cell means measured from the grand mean. 
Holding sigma and "К" constant in the formula for "c", the 
Shrinking tactor is affected py 553. Tne Shrinkixrg factor 
"C" becomes smaller as SSB gets smaller. The predicted 
means are more affected py this condition and the resultant 
wean 1S closer to the grand mean. On tne other hand, as SSB 
increases, "c" increases and the shrinkage is less drastic. 
In effect the James-Stein procedure makes a preliminary 
estimate that all the unobservable means are near the grand 
mean. if tne data supports tne condition wherein tne sample 
üeans are not too far from tne grand mean, then tne esti- 
mates are all snrunk further toward the grand mean. If this 
condition is contradicted, then the Shrinking 15 minimal. 


The distribution of the sample means around the grand mean 


effects tne shrınking factor. The number of means being 
estimated also influences the shrinking factor. Іг "К 
rather large, the Shrinking to the grand mean maybe more 


drastic, [πο ۷۷۶+7٦ 
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TABLE XV 
FORMULAS 
ον] inventory at the beginning of fiscal year t(r 1 τ) 
Let. y number of attritions at any time during year t 
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JAMES-STEIN 


The following steps are utilized to generate James-Stein attrition rates. 


STEP I Use a variance stabalizing transform (Freeman - Tukey). 


De ee еге. 
BL ET (Sin Te т) l) - Sin VETE 1)! 





STEP 2: Form the cell means. the grand mean. SSB, and SSE 


Let. X cell means 
Ee ΛΑ. grand mean 
K 


Ler, number of treatment cells 


STEP 3. Compute the set of James-Stein estimators in the transformed scale 
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STEP 1 Invert the transform to produce the attrition rates "r" 
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Which set of values, 'm' or 'p'", is the better indicator 
of the attrition rate? In order to answer that question one 
would have to know the true attrition rate. This true 
average shall be designated with "TH" іп the main text of 
this paper ana by tne Greek letter theta in the Appendices. 
It is the probability that an officer would leave any time 
and is actually ал unknowable quantity. Although "TH" is 
unooservable, an approximation can be made utilizing 
historical data. 

One method of evaluating the two estimates, i.e., the 
Janes-Stein estimatcr and the observed attrition rate m, is 
by simply counting their closeness to "TH". Consider the 
closer the estimator is to "TH" as a success and the farther 
from "TH" as a failure. The number of successes of one 
estimator could then Ге compared to the number of successes 
of the cther estimator. 

A more quantitative method or comparing the two esti- 
nates is through the total squared error estimation. This 
1s weasured by first determining the actual error of each 
prediction , given by (TH-m) and (Td-p) for each cell. Sach 
of these quantities is then squared and the squared values 
are added up. The observed averages "m" would have a total 
Squared error denoted by Em, whereas the squared error of 
the James-Stein estimators would be denoted by ESSE 
Comparing Em and Ejs by the ratio R=En/Ejs, then the 
James-Stein metnod is "R" times as accurate as the coserved 
attr ı tion rate "ш", 

Employing the ideas borrowed from statistical decision 
theory, the estimators can be compared through a risk func- 
tron: The risk function of a decision rule is the expected 
loss (over the sample) incurred using the decision rule 
۲62.٦93. πε. ο]. The above techniques are similar to the 


squared-error loss method (Ref. 9: p. 9]. 
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Thererore, the risk function is the expected value of 
the squared error for every possible value of "TH". Stein's 
theorem is concerneä with estimation of several unknown 
means (ЕеҒ. 7). No relation between the means need be 
assumed and indeed are assumed to pe independent of one 
another. The historical development of statistical theory 
from Gauss through decision theory argues that the average 
is an admissible estimator as long as there 15 just one or 
two means to be estimated. An estimator for a parameter is 
admissitle if, according to a specified criterion, tnere 
exists nc other estimator that is better than it for all 
possitle values of the parameter. Stein in collaboration 
with Janes showed that when the number of means exceeds two 
that estimating each cf the means by their Own average is an 
inadmissible procedure. No matter what the values of true 
mean, there are estimation rules with smaller total risk 
[Ref. 6]. 

The risk function for the James-Stein estinator is less 
than the risk Luncticn for the sample means irrespective 
what the true values of the means "TH" happen to te. The 
reduction of risk can be substantial, particularly when the 
number of means is relatively iarge. Tne risk of the 
James-Stein estimator is smallest when all the true means 
"TH" are equal. As the true means increase in variation 
from one another the risk of the estimator increases, 
approacaing the value of the observed averages Lut never 
quite equaling it. Ihe James-Stein estimator does supstan- 
tialiy better than the averages only if the true means lie 
near each other. Ihe James-Stein estimator does at least 
marginally better no matter what the true means are (Ref. 8: 
р. 1241. 
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C.  OIHER JAMES-SIEIN MODELS 


The model for the James-Stein estimator used thus far 
shrinks the observed averages to tne grand mean. Ihis is 
not the only possible procedure. Other models ror the esti- 
nator dispense with the grand mean entirely. The observed 
attrition rates do not depend on a choice of origin.  Berore 
Stein discovered his method it was generaliy accepted that 
such "invariant" estimators must be preferable to those 
whose predictions change with each choice of an origin. τς 
the origin or zero is chosen as the grand mean then the 
terms containing the grand mean would be removed from the 
equation. In this case the  James-Stein estimator would be 
р=с (ш). The estimation process is now complicated by the 
fact that the shrinking factor "с" would be different for 
each cell. Tke shrinking factor is dependent on the stan- 
dard deviation of the sample attrition rates. A large stan- 
dard deviation implies a high degree of randomness or 
uncertainty in the sanple data. If the sample attrition is 
large, it can be attributed to random fluctuations rather 
than to an actual large value cf the true mean "TH", Tnus, 
application оғ а smali shrinking factor would reduce the 
value. [Bef. 7: pp. 123-4] 

There is one purpose for which the measured  sauple 
attrition rate may well be superior to the Janes-Stein esti- 
mator: when a single celi is considered in isolation. Ihe 
James-Stein method gives better estimates for a majority of 
Cells and it reduces the total error of estimation for the 
sum of all tne cells. Estimating the true mean for an 
isolated celi by Stein's method creates serious errors when 
that wean has an atypical value. The reduction of the Lisk 
is more predominant in a nomogeneous setting. The inclusion 
of an atypical or ncnhomogeneous mean in the estimation 


process would increase the risk. 


Ss 


de Bayesian versus James-Stein 
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The formula гог the James-Stein is similar to the 
Bayesian equation of Z=m+C(y-0). Here "y" is the sample 
mean, "n" is the mean of the prior distribution. The 
shrinking ractor "C" is different in that it depends on the 
standard deviation of the prior distribution. (Ref. 9: pp. 
99-110] 

The James-Stein procedure, however, has cne impor- 
tant advantage over the pure Bayesian method. The 


James-Stein estimator can be employed without knowledge of 


the pricr distributicn. Tne James-Stein estimator can be 
regarded as an empirical Bayes rule. The empirical Bayes 
approach uses historical data to estimate the prior 


distribution [Bef. 9: p. 117]. 


2. odifications to the James-Stein Process 


Ihe James-Stein estimation process can be modified 
using the Efron-Morris limited translation version of the 
estimator [Ref. 10]. This modification ensures that the 
estimator ος "TH" is not shifted so far from the sample mear 
tnat the estimator is inconsistent with the sample mean. 
The szodification is a compromise between the James-Stein 
estimator and the Maximum Likelihood Estimatcr which has 
good individual properties. The compromise would follow the 
James-Steir rule as closely as possible subject toa fixed 
constraint cr now far the estimator is allowed to deviate 
ρου tne MLE. Such "limited translation estimators" were 
discussed in a Bayesian framework in tne  Efrcn-Morris 
article (Ref. 11). They showed that it is possible to 
consideraply reduce the maximum possible risk for any compo- 
nent while reducing the ensemble savings of the James-Stein 


estimator by as little as 5 or 10 percent [ Ref. 10}. 
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D. VALIDATION TESTS 


A validation test was conducted to evaluate the eifi- 
ciency of tae James-Stein shrinkage estımator. The test was 
conducted as follows: 

1. Select a grade within an Occupational Group to 
examine. The resultant desired data array will be 
three dimensional (years, LOS, MOS). 

2. Let "i" stand for LOS, then 1-0 nnmnnn 

j. Let "3" stand Lor 7 (then values of j‘'s depend on 


which MOS group is being analized). 


4. Let D1j = Incidence Matrix or nonstructural zeroes, 
which is the same for all years. Dij = 1, if cell is 
member or the feasible region; Dij = 0, otherwise. 


Se Let K = number of feasirle cells, 1.e., sum of all 


Dag 
6. Let Yij = number of leavers in cell (i,j). 
7. Let t = 31,...,1; where T = number of years of data 


used to create tne estimator. 
8. Let N1 j = Inventory in cell (173) = 
4ax ((N (t) *N(t* 1)) /2, Y (t)). 
The validation procedures used t-1...4 to compute the empir- 
ical estimates and used t-5,6 for validation purposes. 
Three estimation methods were enpioyed: James-Stein, 


Minimax, and Maximum lixelihcod. 


l. Preliminary Steps: 


NES MENS — m m O F T TD 


The following steps were utilized to prepare the 
data for the validation procedures. These steps included 
the transformation or the given data of leavers and computed 
inventory according to the variance stabilizing arc sine 
transformation listed in Table XV. See Appendix D, Figure 


D.18 for the APL listings of the following functions. 


с0 


Let iSij(t) = Yij(t) BINPREP Nij(t), where "BINPREP" 
prepares the Freeman-Tukey version of the arc sine 
transformation fcr binomial data. 

Iet T = 4, years of data desired. 

Let Jij = Dij JAMES ISij(t), where "JAMES" returns a 
James-Stein estimator for the means of the cells in the 
last two dimensicns ot "IS" while being screened by the 
incidence matrix "D" of nonstructural zeroes. 

Let Pij = Jij BINCONV Nij(t), where "BINCONV" inverts 
the arc Sine transformation used in preliminary step 
cne. 


The rollowing procedures were utilized to validate 


the effectiveness of the the James-Stein estimation process. 


Бес” 15491) = 1Sij(I+1), wnich is assumed со be 
distributed Normal(Jij,unknown variance). 

let NOR = І5%137 - Jij, where Jij was derived in 
preliminary step three. 

Use the distribution fitting capapility of GRAFSTAT or 
any other comparable graphical display package, to 
compare NOR to the Normal distribution. Desired output 
is a comparison histogram of the data to a fitted 
Normal density, compare sample distribution to ritted 
Normal distribution, normal probability  Q-Q plot with 
data, and possibly the survivability curve fit. These 
plots can all be plotted on one Sneet of output. 
Compute test Statistics comparing data fitted to a 
Normal distribution. 

Compute the figure of merit from the test statistics 
using the sum of the square of the Sampie mean and the 


square of the Sagpfle standard deviation. 
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The following procedures were utilized tc compare 
the James-Stein estimator to the Minimax estimator. See 
Table XV for the formula for the Minimax eStimatcr. An 
simple example of the possible values resulting from the use 
of this formula is provided in Table XVI, to illustrate the 
shrinkage characteristics inherent in its use for the 
purpose of projecting rates which have extreme empirical 


values. 
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TABLE XVI 
MINIMAX EXAMPLE VALUES 


MINIMA 

LEAVERS INVENTORY ESTIMA 
0 0 . 50 

1 25 

1 5 
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See Appendix D, Figure D.19 for tne APL listings of tnis 
Funetıon. 

е Tet PMM3JJ = 1η MINMAX Nij(t), where "MINYAX" 
returrs the MINMAX estimates ror the binonial. 

e Let ISM = Nij ARCSIN PMMij, where "ARCSIN" returns the 
inverse sine transformation for use when the success 
probability PMMij is estimated directly  (e.g., by 
I NMAX). 

e let NORMM = IS* - ISM, where IS* is from validation 


procedures step cne. 
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e Use the distribution fitting capability of GRAFSTAT or 
any other comparable graphical display package, to 
compare NORMM to the Normal distribution similar to the 
validation procedure three. 

e Compute test statistics comparing даса fitted to a 
Normal distribution. 

e Compute the figure of merit from the test statistics 
using the sum Of the square of the sample mean and the 


square of the sample standard deviation. 
4. Maximum Likelihood 


Ihe natural empirical estimator without shrinkage is 
the Maximum Likelihood estimator. Comparison of the 
James-Stein process to the Maximum Likelihood process will 


illustrate the efficiency of tne James-Stein estimation 


process. 
e Iet Y = sum of Yij(t) for t-1 to T. Where T = 4. 
e let N = sum of Nij(t) for t=1 to T. Where T = 4. 


e jet ISAij = Y BINPREP N. 

е let NORA = IS+ = ISAilj, where IS* is from the valida- 
ticn procedures step one. 

e Use the distribution fitting capaoility of GRAFSTAT or 
any otner comparable graphical display package, to 
compare NORA to the Normal axzstribution similar to the 
validation procedure three. 

e Compute test statisticS comparing data fitted to a 
Normal distribution. 

e Compute the figure of merit from the test statistics 
using the sum of tne square of the sample mean and tne 


Square of the sample standard deviation. 
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E. DATA ANALYSIS 


A data analysis was conducted utilizing the three 
comparison mentioned above. 


— | DA. F ee ee ч- ч s s cu 


Four years of data (1977-1930) was used to construct 


the James-Stein estimator. The leavers and inventory were 
transformed with the Freeman-Tukey transformation. The 
projected year's (1981-1983) leavers and inventory were 


constructed and also transformed vita the  Freeman-Tukey 
transformation. The difference between the  James-Stein 
estimator and the projected years data was plotted to 
compare tne rit to a Normal distribution. The figure of 
nerit (РОМ) was computed from the inputs taken from the 


Statistical table. 
FOM= (Squared mean) * (variance). 


Table XVII illustrates the results of the computa- 
tion cf the James-Stein, Minimax, and Maximum Likelihood 
BOM SS The data organization format of the table 45 
provided for ease of interpretation of the table. The table 
demonstates the calculation of the FOM for the years 1981, 
1982, and 1983. The separate cases evaluated are identified 
by rowS which represent the grades of First Lieutenant and 
Lieutenant Colonel in each of the four MOS groups utilized 


ІП “Ете project: 


2. Zhe Minimax Estimator 


The leavers and inventory of tae four years 
(1977-1980) were used to construct a minimax estimator. See 
Table XV for the formula utilized to derive this estimator. 
The result was compared to the fit of the Normal distribu- 


tion іп а similar manner as discussed above. The FOM was 
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computed as shown abcve and further illustrated in Table 
XVII. 


3. Maximum Likelihood 


The last comparison was using the leavers and inven- 
tory of the four years (1977-1980) and only using the 
Freeman-Tukey transformation. The result is the Maximum 
Likelehood estimator in the transformed space. The algo- 
rithm is the same as James-Stein estimator described above, 
only the James-Stein calculation is by-passed. The 
projected years data was used to compare the resultant 
ditference with the fit to the Normal. The FOM was computed 


as shcwn above and further illustrated in Table XVII. 


4. Figure of Merit 


The figure of merit is representative of the effi- 
ciency of the estimator. Since the FOM is computed in the 
transformed and Normalized scale, the smaller the FOM, the 
petter the estimate. The ratios illustrated in columns 3 
and 5 of Table XVII illustrate (in most cases) how many 
times larger the other estimating technique's FOM is when 
compared to the James-Stein estimator technique. One can 
note that in one case the use of the James-Stein estimator 
technique resulted in an improvement of as much as 30 times 
over the Maximum Likelihood process and 22 times over the 
Minimax process. These results are typical or other cases 


executed by tne author, but not snown for redundancy sake. 


F. ADDITIONAL COMPARISONS 


Ihe above risk ccmparisons were all performed in the 
transformed (Freeman-Tukey) space. Since this distorts the 
scale it seems wise to perform some additional comparisons 


(sun o£ squares error) іп the original (untranstorned) 
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TABLE XVII 
FIGURE OF MERIT COMPARISONS 


*#OHSHEEXEKEK DATA ORGANIZATION FORMAT жжжжжжжжжжж 
*o** COLUMN  x*&x* *** REMARK  *** 
JAMES FOM of data using James-Stein 
MAXLIKE FOM or Maximum Likelinood 
MAXLIKE/JAMES Ratio of Maxlike divided by James 
Mi NIMAX POM of data using Minimax 
MINIMAX/JAMES Ratio of Minimax divided by James 
*** ROW жжж *** REMARK **% 
1 IST LT UB Aviation ۱ 
2 1ST LT in Combat Service Support 
3 1ST LT in Combat Supe E 
5 IST ET 8ص1‎ 677٤ 
S LICOL ın Aviation 
6 LTCOL in Combat Service Support 
7 LTCOL in Combat Support 
8 LTCOL in Ground Combat 
ΚΑΚᾺ THESE ARE THE FIGURES OF MERIT FOR 1981 Жа 
MAXLIKE MINIMAX 
RCW JAMES MAXLIKE JAMES MINIMAX JAMES 
1 1.69 51.47 30.46 3272668 2229. 
2 0.71 1. οὐ O 0.75 1.10 
3 1.20 21019 32.03 2.10 1.78 
4 1. 19 7د‎ q. 59 dad d 320 
Б 3.05 13.64 4.47 В 2:92 
6 (7522 OS UE 1.91 ο 3 1205 
7 OSM 2.34 7.55 πο 3.74 
8 0279 6.81 0.02 22768 4.72 
¥ kk ¥ THESE ARE THE FIGURES OF MERIT FOR 1982 7 
MAXLIKE/ MINIMAX 
RCH JAMES MAXLIKE JAMES MINIMAX JAMES 
1 4.37 47.05 19272 30. 29 7485 
2 0-79 1.55 07 0.76 1.01 
3 1.53 3. /1 2. 42 20۲ 1.43 
y . 5.84 Sis ors 4.21 2.28 
Е 380 17225 238 192.08 2-26 
6 052 3 0.54 2ً ο 02093 1.43 
7 O. 59 2250 4.41 1.33 2.42 
8 JS 8. 34 T a 4.92 435 
жажж THESE ARE IHE FIGURES OF MERIT ۶ 0۶۰ ۱۹۹33 52 
MAXLIKE/ MINIMAX 
ROW JAMES ΛΕΣ JAMES MINIMAX JAMES 
1 4.04 43.16 η UNOS σ ο 
2 0.91 1.64 1-80 2:306 05 
3 1.17 3. 67 3.14 2.02 1.73 
4 285 5.67 1. 3 4.21 1.48 
3 5.49 11. 93 2:217 8.44 1. 94 
€ 0.20 0. 56 23855 0.34 . 
7 0.62 2.76 4.45 1.55 2.47 
8 1257 8253 5.45 5.16 2524 
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scales. Accordingly the deviations of the empirical attri- 


tion rates in the validation time frame from the estimates 


(both James-Stein and Maximum Likelihood) are computed and 


compared. Specifics follow: 


1. Subtest 


The preliminary steps are as follows: 


Actual rates: empirical attrition rates are the central 
attrition rates fcr the prcjected period. 

James-Stein Projected rates: projected attrition rates 
utilizing the James-Steln estimation process previously 
defined. 

let, MLE= Maximum Likelihood Estimates, the historical 
empirical attrition rates are assumed to remain the 
same throughout the projected period. 

Tet, SEP- The squared differences  Letween the 
James-Stein projected attrition rates and the actual 
Ia tes. 

let, SEM= The squared differences between the MLE and 
the actual rates. 

let, SSE(MLE)= Ihe sum of squared differences of SEM. 
let, SSE(P)= The sum of squared differences of SEP. 
Let, SEP<SEM be the number of cells for which SEP was 
less tnan SEM. (For convenience of column heading, 
inequalities are allowed to represent numbers.) 

let, SEM<SEP be the number of cells ror which SEA was 
less than SEP. 


KOS Groups evaluated were: CS - Combat Support Group; 
CC - Ground Combat 
banks evaluated кеге: ІТ - First Lieutenant; LICOL - 


Iieutenant Colonel 


| 


First it was assumed that all cells in tne feasible 


region were eligible fcr forecasting tne projected attrition 


rates. The results ог subtest 1 are illustrated in Table 


o 7 


7۷.۰۲ The sum of sguared differences for the MLE and the 
James-Stein projected rates fcr the first and last groups 
are competitive. One has to keep in mind that the 
Jame-Stein estimation process will result in a value greater 
than zero given a cell with an original value of zero. 
Therefore, wnen ail tne cells in the feasible region are 
considered eligible fcr projection, the empirical values of 
zero result in a value closer to the grand mean. The value 
or the sums of the differences between the James-Stein esti- 
mation process and the actual rates can be understandably 
larger than the sums of the differences between the MLz and 
the actual rates when there are numerous zero cells 


originally in the feasible region. 
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TABLE XVIII 
SUBIEST 1 


(Қа کک‎ p 


Feasirle 
Grcup Rank Cells SEP<SEM SEM<SEP SSE(MLE) 55Е(р) 
CS LT 57 15 ША 52 1.91 
GC τι 57 22 35 Ú. 33 3.99 
CS II OL 48 12 36 1.29 2.39 
GC ТЕСЕ 48 21 21 2.05 0.66 
LIII —— — — іл созы ты сы са eee 


2. Subtest 2 


A second comparison assumed a conditioning that 
considered only the celiS in which there was an actual 
attrition rate greater than zero. In other words, what was 
the projected attriticn rate given an actual  lcss occured? 


This created a subfeasiole region whicn was smaller than the 


whole feasible region. The results of subtest 2 are illus- 
trated in Taole XIX. The sum cf squared differences for tne 
James-Stein projected rates in this subtest are much lower 
than the MLE. 


| 
| TABLE XIX | 
| 
| SUBTEST 2 | 
| suD- _ | 
| Feasicle | 
| Grcup Rank Cells ЗЕБСЗЕП 5ЕМ<5ЕР SSE(MLE) SSE (P) | 
CE LT 26 15 11 1:52 0227 
GC γα 29 22 7 0233 0.10 | 
es L TICO 16 12 4 1229 0.14 | 
| GC ΘΙ 28 21 7 2905 ΙΖ | 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Pas EEE Tr u СЫЗ کیوے-سوو‎ TE. ON UNA жә» «и» ақш» um Qan سے‎ CA Da s ee ee ce as s Gm 


A. SUMMARY 


The purpose of this study was to demonstrate the appli- 
cation of the James-Stein and other shrinkage type parameter 
estimaticn schemes for the goal of generating Manpower loss 
rates within the USMC officer force structure. This thesis 
contains comparisons of performance of James-Stein, Mininax, 
and Maximum Likelihocd estimations of Marine Corps officer 
attrition rates for select cell aggregates. The very large 
number of cells within the USMC officer force structure 
leads tc the conditicn that empirical attrition rates are 
unstable. This problem is compounded by the fact that the 
cell protabilities are small. Most rates are less than 10% 
and virtualıy none are greater than 20% at disaggregated 
levels. Further difficulties are preseat because some oL 
the inventory cells are empty for structural reason while 
otners are empty by chance. Theretore, the small inventory 
cells draw especial attention. It has been illustrated 
within tris project [παῖ improvement can be attained by 
application of the James-Stein and Minimax Shrinkage methods 
ratner than the more natural Maximum Likelihocd estimation 
process. It is impcrtant to note tnat the James-Stein and 
Minimax Snrinkage schemes seem to compete for appropriate 


attrition rate generation. 


В.  CCNCIUSIONS 


The shrinkage schemes emplcyed herein offer powerful and 
useful methods for generating attrition rates. This employ- 
ment should lead to much lower costs that accrue from errors 
in manpower planning. However, the particular way to 


exploit them is yet tc be determined. 
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The application of the James-Stein shrinkage technique 
has exposed a number of problem areas which could be 
expanded. 

The aggregation problen. Aggregation of cells (cross 
classified Ly grade, Military Occupational Specialty, length 
of service, etc) intc sets whose attrition behavior is hono- 
geneous and of low internal variability should lead to the 
most useful overall performance. One method, not attempted 
by the author, is to aggregate over MOS types in which the 
grade structure is similar. For example, some  MCS's are 
restricted to warrant officers only, warrant officers and 
limited duty officers only, unrestricted regulary/reserve 
ofticers only, etc. Another method would be to aggregate 
over certain length of Service periods. For example, 1-4 
years of service, 5-10 years of service, 10-20 years of 
service, ard 21-30 years of service. Thus examining the 
attriticn rates in relationship to preakpoints and lcngevity 
in the officer's career pattern. 

The limited translation problem. The shrinkage estima- 
tors translate tne rates toward the agyregate mean. Such 
translations can appear to be excessive in instances of 
extreme rates. This suggests that techniques that limit the 
translation could be profitably applied. Current literature 
contains limited translation methods wnich could be 
examined. 

The yearly update problen. The estimated attrıtion 
rates are based upon data from several recent years. AS 
each year produces a new set orf experiences there is need 
for a policy to include new data in the estimation scheme 
and phase out the influence of the distant past. 

The validation problen. The smallness of the rates 
(1.e., 0-20%) places the current project inr a range which 
nas nct yet been treated successfully in the literature. 
For this reason it is especially important to apply cross 


validation procedures to score the efficacy of the metnods. 
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The multinomial estimaticn problen. The number of 
attritions in a cell are placed into a number of (currently 
eight) disjoint categories. This produces a number of small 
sample size multinomial probability estimation problems. 
Metnodologies for managing such problems have attracted 
attenticn in the recent literature. Their usefulness could 
be further examined in regard to the officer attrition rate 


estimaticn process 


С. RECOMMENDATIONS 


It is recommended that further studies pursue the prob- 
lens identified above. It is suggested that a new data tape 
be provided to the Naval Postgraduate School which allows 
the distinction between restricted and unrestricted  orfi- 


cers. Additionally, the data for fiscal years 1984 and 1985 


needs to be provided in the near future. In order to 
provide coLtinued service, a regularly scheduled yearly 
updated data set should be torwarded to the Naval 


Postgraduate Scnool when available. 


CENTRAL ATTRITION RATE GRAPHICS 


This appendix contains graphical illustrations oí the 
central attrition rate for USMC officers. 

Figures A.1 to A.10 contain the central attrition rates 
for each Military Occupational Specialty distributed over 
the length of service period from 0 to 30. These figures 
illustrate the aggregate central attrition rate for all 
eight types of losses over the six years trom 1977 to 1982. 
These figures demonstrate the  breaxpoints within each MOS 
over the distribution of length of service. 

Figures A.11 to A.20 contain the central attrition rates 
ror each Military Occupational Speciality distributed over 
the grades trom Warrant Officer to Colonel (latelled 0...9 
on the akscissa). These figures aiso illustrate the aggre- 
gate central attrition rate ror all eignt types of losses 
over the same six years. 

Figures A.21 to A.28 contain tae central attrition rates 
for each Military Occupational Specialty distributed over 
the years from 1977 to 1982. Each figure has rive MOS 
graphs and are distinguishacle by the type of line and 
symbol pcints for each separate 405. These figures illus- 
trate the aggregate central attrition rate fcr all  eignt 
types ot losses over the six years from 1977 to 1982. 

Figure A.29 contains the central attrition rates ror the 
Aviation Group for the grades from Second Lieutenant through 
Lieutenant Colonel for tae years 1977 to 1982. Figure A.30 
contains the central attrition rates for the Combat Support 
Group for the grades from Second Lieutenant to Lieutenant 
Colonel ror tne years 1977 to 1952. Figure A.31 contains 


the central attrition rates for the Ground Comtat Group for 
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the grades from Second Lieutenant to Lieutenant Colonel for 
the years 1977 to 1982. Figure &TOTMOS contains the central 
attrition rates for the Total MOS's for the grades from 
Second Lieutenant to Lieutenant Colonel for the years 1977 
το VIEL 
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Figure A.29 Central Attrition Rates for the Aviation Group. 
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Figure A.31 Central Attrition Rates 
for the Ground Combat Group. 


105 


Lomo سے ےہ سے جو‎ Bien m SE a nae SE CS a A eter PA ыл}. (Шы илии». Шин». A | 





APPENDIX B 


JABES- STEIN ESTIMATOR ALGORITHM 


This appendix contains the detailed algorithm for the 
Janes-Stein estimator of attrition rates for the USMC 
officer force structure. The information provided includes: 

e Assumptions. 
e Notation. 


ə Steps to generate the James-Stein attrition rates. 
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STEP 4: Compute the set of James-Stein estimators in the transformed scale. 
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Meee Le STOCHASTIC MODEL FOR CENTRAL DATA 


This section contains the nomenclature and characteristics of a simple stochastic model for 
central data. It deals with the joint distribution of cell inventory and number of attritions. 


The following nomenclature will be defined in this section as follows: 
N = Inventory during period, which is distributed Poisson ( A ) 


Y = Number of losses. 


m = Central attrition rate = — 


> 
The conditional probability of losses given inventory is: 
p(¥iNen)= (7) evti- e)o7* Un 


The Joint probability function is: 
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D SUPPORTING DISTRIBUTIONS. 


The following sections of this appendix contain a series of proofs to the following statements 
based on the assumption: Y is distributed Normal (©,/) . where I is the Identity matrix. 


SUI 


j| ہہ‎ is distributed Normal (| © b 

nee = Y - oe = Y- = is distributed Normal (0,7 жылы. 
ОР x ERE rer]? 

See section C for the following statements: 

ii) V=/¡Z] | % is distributed xx -1) - 

ο =| UV ετσι". 


PROOFS. 
1) U is independent of all components of Z. In other words all Cov(U,Z)= 0. Proof is 


as follows. 


Cos (U ,Z ) = ہتشر وج‎ = Cov (U, Y )-—— Cov (UU) 
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C DISTRIBUTION " V" 


This section is an continuation of section B and shows V - Z Z = 
l 


buted Y(k -1) - у) Normal [0.1)°. 
l 


D а 
27 Se distri- 


PROOFS. 
iii) First, construct an orthonormal transformation P, where P is K by K dimensional and 
P’P = I, where I is the identity matrix. 
Let px be the last row of P and be defined by 
7 Θ 
770. Θ 
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` PK = 1 
Note that this can be done because an orthonormal tranformation can always be completed. 


Now let Y = PZ and calculate 
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Therefore X4 = O random variable. and the space is K-1 vice K dimensional. 


Next calculate 
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D. 


NON-CENTRAL CHI SQUARE 


This section considers the Non-Central Chi Square with n degrees of freedom. 


Letting 


ieee, Y, be tid Normal (0,1). then x4,= У; Y * with density: 
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But if the means are not all zero then the result is the Non-Central Chi Square random vari- 
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of V ıs as follows. 
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With n=1 and letting Y be distributed Normal (0,1) and 







































































V = Y”, then the density 
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E. DISTRIBUTION of W and EXPECTATION of W^! 


This section illustrates the characteristics of the distribution "W". which ts the convolution 
of f ,.2(8 ) (non-central chi square) and the density of Vk Eon uis noted that W is a 
non-central chi square random variable with K degrees of freedom and non-centrality parameter 
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F. LOSS AND RISK FUNCTIONS 
This section evaluates the Loss and Risk functions. Considering the column vector Y. which 


is distributed Normal (©./), these functions are as follows. 


LOSS - L(a.0) - (a-O)(a-O)- а-ө 2 


Risk is expected loss. and for the maximum likelihood estimator, a — Y and 


RISK - E; L(Y,.O) | - E(Y- oy (Y-O)- E(, Y-O0 i : 3) = K 


4 


5 
Consider another function of Y. Use Ie ууз Y Y,* which is distributed 
ڈور‎ 
Non Central Ur) See section D for the properties of the Non Central Chi Square Distribution. 
The James-Stein estimator has structure: 


b 


a (Y) = (| ——— } ¥ 
ο ο.» —3 ) 
It's Risk function is derived as follows: 
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] E τ᾽ 
where U and W are defined in section B of this appendix. 
See section G in this appendix for an evaluation of El Using that result now leads to 


the following: 
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ESEVALUATION OF )اط‎ - ) 


TEN U 
The goal of this section is to evaluate E| ЕТЕ ho 


First, find the joint distribution of U and W. 


fuvlu,v) = fulu)fylv) 
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ο ΕΥ VERSES YEARLY RATES 


0۵ 

This section provides a recommended method to convert yearly attrition rates to monthly 
rates. 

Let r = yearly attrition rates. 

Let r,, = monthly attrition rates. 


Assume independence from month to month. Thus 12 consecutive survivor months implies 


a Survivor year. 
12 
l="=(1-r,) 


1 


кугу! ле zer) 
ake 
η 


Tm 


Thus. the above equation relates the vearly attrition rates to the monthly attrition rates 


ιο. 





MSSCOUN JCL A 
//COUNTS JOB (2209,5555),'MAJ D. D. TUCKER',CLASS=B 
// *MAIN ORG zNPGVM1.2209P 
//* t*sustttsueusanseseses«e NOTE TO СЕКС +++» жя я жактан 
//* THIS WILL MOVE THE DATA SET COUNTS FROM MASS STORAGE TO THE 
//* МУ5004 DISK. FIRST CHANGE THE XXXX'S TO YOUR USER ID. 
//* CHANGE YYYY TO YOUR PROJECT NUMBER, THEN SUBMIT THE PROGRAM 
//* TO MVS. WHEN THE SYSTEM RETURNS A COMPLETION MESSAGE, USE 
//* THE GETMVS SYSTEM EXEC TO COPY THE DATA SET TO YOUR A-DISK 
//ж AF YOU WISH TO ACCESS IT FROM CMS. BE ADVISED THAT THIS 
//* FILE IS HUGE AND WILL TAKE UP NEARLY 4 CYLINOERS. 
//* 
//ж USE DSN= $2209 COUNTS IN GETMVS 
//* دبع راچ عو جج عو رجہ رھ تھے‎ RARER RE REE Καὶ ھچ جو و اد وس‎ αὶ Καὶ ES 
// EXEC ΡΟΜΞΙΕΒΟΕΝΕΕ 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD DUMMY 
//SYSUT1 DOD DSN=MSS.F0968.COUNTS,DISP=OLD 
Poise sO DsN=Ss2209 COUNTS, VOL=SER=MVSO04,UNIT=3350, 
LA SPACE=(CYL,(2,2),RLSE),DISP=(NEW,KEEP) 
/* 


ws A ----- ge (QQ ست .ری‎ en اا‎ ee 
L eee ہے‎ eee eee ee Oe ae eee nm nn nn nn о ου μμ ουν a, CEE gal 
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Figure D.1 MSSCOUN JCL Al. 


B. |. DATA MANIPULATION PROGRAMS 


A technique utilized zy tne author to reduce the depen- 
dency of manipulating such a large data rile, involved 
Creating a rile of indices corresponding to pertinent Jata 
in the original summary data file. A Set or indices was 
Created for the inventories and losses identified ror each 


fiscal year. 


The set indices for tne inventory cat ре created 

usirg tre follow1ry procedures. 
e Separate the sugmary data file by fiscal year data and 
Create a new file tor each riscal year. For example, 
۷۰۰۰٠۱٠٠٠٠٠۰٠۰ ۸٣۲۹۰ ت‎ Tor” fiscal Year 1977 data on tne 


teuporary C-disk. 





x 
۱ 


TABLE XX 
RECORDS PER FISCAL YEAR 


Fiscal Year Number oc kecords 
1977 2203 
1978 2191 
1979 2337 
1980 290 
1981 ΙΙ 
1982 2324 
1983 2330 


cm Se Se سس‎ Re MS RD ee ہن کہہے‎ 


eo ee er μασ IA as == 


FILE: NUMBER EXEC Al 


TRACE OFF 
x TO EXECUTE THIS EXEC TYPE "NUMBER <YEAR>" 
x NHERE YEAR IS OME OF 77-83. 


GLOBAL TXTLIB VALTLIB VFORTLIB IMSLSP NONIMSL CMSLIB 
GLOBAL LOADLIB VFLODLIB 

FORTVS INV CLVLO77)) 

* FORTVS 81 (LVL(77)) 

CE TERM LINESIZE T55 

SA = SCONCAT OF COUNT &1 

&B = &CONCAT OF INV 41 

FILEDEF 08 DISK &A DATA C (RECFM F LRECL 53 PERM 
ΕΠΕ ΠΚ 3B ARRAY C CRECFM F LRECL 9 PERM 
&TRACE ON 

LOAD INV (START 

X LOAD 81 (START 

ERASE INV LISTING 

ERASE INV TEXT 

ERASE LOAD MAP 

&EXIT 


— = AO A A e шшш ہے کت‎ =w— = 


Figure D.3 NUMBER EXEC A1. 
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«Ine CMS file "INVXX ARRAY C", саг be read into an AP L 
workspace and assigned to a numeric variable by 
executing the user friendly "CMSAPL EXEC AU. See 
Figures D.5 and D.6. Note that "CMSAPL" has a provi- 
sicn to add a variable to an existing APL workspace or 
to create a new workspace each time  "CMSAFL" is 
executed. 

The file of indices created in this manner can be 
utilized to generate matrices in the APL workspace. The two 
dimensional array for each fiscal year's inventory is char- 
acterized by a length subject to the number of records in 
the summary data file ror a particular fiscal year in waich 
there is a nonzero inventory count. The columns of the 


records in the index file is described 1n Table XXI 


2-2 loss indices 


Ine set indices гог the losses can be created іп а 
manner similar to the indices for the inventory by using tne 
following procedures. 

e Separate the summary data file ру fiscal year data and 
create a new file for each fiscal year. For exanple, 
"705577 DATA C" for fiscal year 19277 data on the tempo- 
гагу C-disk. 

«Ех еспе спе “LOSS EXEC A1" cy typing “LOSS <Xx>" where 
XX 15 the desired fiscai year (such as 77 for fiscal 
year 1977 data). See Figure р.7. This executive will 
ο l ου Еос ап program» “LOSSES FORTRAN Al" which 
reads the data trom the file "COUNT77 DATA С". See 
Figure D.8. 

e The output from “LOSSES FORTRAN A1" wili be sent to a 
СМ5 file  "LOSSXX ARRAY C", wnere XX is the desired 
fiscal year (same as above). 

e The CMS tile "LCSSXX АККАҮ С", сап be read into an APL 
workspace using similar methods as tne inventory tile 


described above. 
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m 
| 


EP a a m es = m сш 


&STACK N 

&STACK )WSID MYSPACE 
&STACK )SAVE 

&STACK )OFF HOLD 

EXEC APL 

-SKIP «ТҮРЕ 29222 IS THIS A NEW WORKSPACE? Y OR N 
XREAD VAR &NEW 

&IF &NEW EQ Y &GOTO -JUMP 
CP TERMINAL APL ON 

&STACK )LOAD &WKS 

&STACK )PCOPY MYSPACE &A 
&STACK )SAVE 

XSTACK )DROP MYSPACE 
&STACK )OFF HOLD 

EXEC APL 

&EXIT 98 

-JUMP CP TERMINAL APL ON 
«STACK )WSID &WKS 

&STACK )PCOPY MYSPACE XA 


&STACK )SAVE 
«STACK )DROP "MYSPACE 


&STACK JOFF HOLD 

EXEC APL 

& EXIT 99 

—TELL &TYPE YOU HAVE ENTERED TOO MANY OR NOT ENOUGH ENTRIES 
ТҮРЕ OR THE FILE THAT YOU WANT TO BE READ INTO APL. 

ἘΤΥΡῈ YOU NEED TO BEGIN AGAIN 

«TYPE 

TYPE ENTER: CMSAPL 

XEXIT 100 

-ERROR &TYPE &1 &2 &3 DOES NOT EXIS ON YOUR XCONCAT OF &3 -DISK 
&TYPE CHECK YOUR FLIST AND THEN BEGIN AGAIN 

&TYPE BY ENTERING CMSAPL 

«EXIT 101 


| 
| 
| 
۱ 
۱ 
| 
| 
| 
| 
| 
| 


Figure D.6 CMSAPL EXEC А1 (CONT.). 
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PROGRAM LOSSES 


eT RS ee Tee eT‏ رود کے ے کڈ ٹک 


ο ο ο ο ϱΏϱΏϱϱϱ NESCREATES THE INDEX ARRAY FOR COSS DATA ELEMENT 


SUBROUTINE SUM (INDEX, I,J,K, L, M, NUM) 
INTEGER IND LM, NUM 
DEK IF 3: 100666 00} + (K* 1000000) + (L*10000 ) * (M* LOOO ) +IUM 
WRITE (9,500 
500 FORMAT (110) 
RETURN 


ANA 


INDE 


A 
b 


. 1 
| 
| : 
C *** REFORMATS DATA FROM 'COUNTS' FILE INTO À TWO ٭٭٭>‎ 
С жж DIMENTIONAL ARRAY OF LOSSES FOR EACH YEAR. ыды 
С ... HE OUTPUT IS A SET OF INDICES FOR MATRIX kex 
ς τ. FORMATION IN ABS E 
| : τ T دک‎ INDEX Eo | 
| С .πα 1 1 YEAR 5, πλ zex 
С kee 2,3 2 MOS es 9) -xx { 
LE ee к NT 2 | 
O ΕΞ (O Ll... 
С з 8,9,10 3 NUMB В” бё eo cs | 
σ хжж USE "LOSS <YEAR>" EXECUTIVE xx 
C κ ἡ κα c dg ode ose e ose x πε πε x Ww Ww Ж w x Ww wv w جو‎ IF KKK KK 
| E INTEGER GEE MOS, YEAR, LOSS, LOS, IYR 
| С THE DO LOOPS CORRESPOND TO THE NUMBER ОЕ RECORDS FOR YEARS 77-83 | 
| c THE FIRST ONE IS THE FULL DATA SET, THE COND IS FOR TESTING | 
Ç DO 200 [=1,16093 
С 96 200-1 50 
| 2 2396-200121 3 
E Dnou200 11.2231 
E D68200: S 1° 2337 
| С DO 200 [=1,2351 
E 290-2003121, 7 
C Do02008121 2124 
Е BO 200 11230 
ς ** L1 - L8 IS LOSS TYPE CODE 1 - 8 
101 READ(8, ο. NMOSNGRASDE τος NV L1,L2,L3,L4,L5/L6,L7,L8 
poe FORMAT (412,915) 
| ς ο τα ο στο THE INDEX  1=77,2=78,...,7=83 | 
| ΤΥ YES 76 
E moss = 109000060 * [YR 
' EEL TC O EAL SUM (IN, LOSS, “MOS, ,GRADE,LOS,1,L1 | 
| IF(L2.6T.0) CALL SUM/ IN, LOSS.MOS,GRADE,LoS.2,L2 
| ШЕ(ІЗ ӘЗ ИЕСІН ο 1Ν τοςς Μος GRADE. LOS, 3,L3 | 
| herbs επ σολ SUOMI IN, LOSS, MOS, GRADE, LOS, 3, L4 
REE SC τρ σα LOSS,MOS,CRADE, 5,5, [9 
ПЕ πο GE oO O ο IN LOSS, MOS CRADE.L08,0, Lo 
ENL? ОТОО του a le BDE егы 
Ιεί58.ΟΣ 0) CALL SUML EN. LOSS, MOS, GRADE, LOS, 3,L3 | 
| 200 CONTINUE 
| END | 
| | 
| l 
| اأ‎ 


Figure D.8 LOSSES FORTRAN Al. 


The file of indices tor tne losses is similar to the 
Inventory tile. The two dimensional array ror each riscal 
year's loss is characterizea by a length subject to tne 
number cf records in the Suamary data rile tor a particular 
rlscal year ın Whicu there 1S a nonzero loss count ror each 


Loss туре. It should be Lotel that a record in the original 
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a. Function GETINV 


The function to create the inventory Latrix tor 
each desired year is "GETINV". ее ص۶‎ ٰ 0 "GETINV" 
utilizes the function "INVNATX INVXÉ",  whicn interprets the 
inventory index matrix “INVXX" and creates tke inventory 
matrix for riscal year "KX". See Figure D. 10. The resul- 
КОПТ Matrix is "IXX" for fiscal year "xx". The function 
"INVMATX" could create a matrix of the tollowing dimension 
7x4Ux10x31 ror 7 years, 40 MOS's, 10 grades, and 31 LOS's. 
However, due to limited workspace, the digension of 40x31x10 


for 40 MOS's, 31 LOS's, and 10 grades was commonly utilized. 





ie‏ د I eR Ie SE m‏ ي ڪڪ ی 


V GETINV 
[1] a GET THE NUMBER OF YEARS OF INVENTORY DESIRED. 
[2] a COMMENT OUT YEARS NOT USED IMMEDIATELY. 

[3] 177—INVMATX INV77 

[4] 178-INVMATX INV78 

[5]  179-INVMATX INV79 

[6] I80-INVMATX INVSO 

[7] 181 -INVMATX INV81 

[8] a 182-INVMATX INV82 

[9] a 183-INVMATX INV83 

[10] ' SHAPE OF 177 IS ' 

υπ 


Kan rn nn a ie en m. i ے ے‎ O asan لے سے‎ 


Figure D.9 APL Function GETINV. 


L. Group Inventory Functions 


There is a separate function to take a combined 
matrix for each or the MOS groups. fmeted in Table XXIII 
are the JToups, the function used to create the desired 
Matrix, and the shape of the resultant matrix. Figure D.11 


is a representation of the programing technigue used ror the 





— a Q raw a m T onm rss s sa. cd a mr 


V GETGC;K 


κ 
| [1] в THIS GETS THE INVENTORY MATRIX FOR THE GROUND 
| [2] a COMBAT OCCUPATION GROUP. 
| [3] a SHAPE IS (MOS;YEARS;LOS;GRADE) 
| [4] a «e. CHECK THE 2ND INDEX FOR THE CORRECT 
[5] a «+. NUMBER OF YEARS DESIRED. 
| [6] GC-(3 4 70 
| [7] a THIS GETS THE MOS: INFANTRY 
| [3] GC[1:::]-GETMOS 3 
| [9] a THIS GETS THE MOS: ARTILLERY 
[10] GC[2:::]-GETMOS 5 
| [11] a THIS GETS THE MOS: TANKS AND AMPHIB 
| [12] GC[3;;;]=GETMOS 10 
o 
Figure D.11 APL Function GETGC. 
A 
| V Z-GETMOS X;J;Y 
| [1] a GET THE CENTRAL DATA NUMBER OF INVENTORY 
| [2] a FOR THE YEARS DESIRED FOR A PARTICULAR MOS. 
| [3] a X= MOS DESIRED. CHANGE THE INDEX FOR Z IF 
| [4] a YEARS ARE NOT SEQUENTIAL. 
| [5] Jexti 
[6] ا‎ 31 10)00 
| (7) 2Z[1:)(Y-(77[3;,]+178[3;:))=2351.77[3;:D)) 
i ο] ο ο πο πο 2) B :))) 
| [9] 2[3;,]J-((ve(179[3;:]+180[3;,)) =2)M(179[3;;]») 
[10] z[#: ee. (180[4,:]+181[3;:])= > (So p) 
| (11) az[5::]-(Ye(d81(;:]+18203::) 218 1(4::))) 
1 ШЕШЕ Пар ELL DD OCOD 


o == 2 


Figure D.12 APL Function GETMOS. 


To create the inventory matrix used Lor 
James-Stein runction the roilowiay functions were used: 


e GETNG - creates the "NG' matrix ror orouna Comtat. 
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= 


the 


See 


Figure D.13 for an example of the program utilized to 


create the the Ground Combat inventory matrix. 


The 





a. ۲۰56٦190 СЕТІОЅ5 


The functicn to create the loss matrix гог each 
desired year is "GETLOSS. See Figure 0.14. MGPDBDSST 
utilizes the function "MATRIX LOSSXX",  wLicn interprets the 
loss index matrix "LCSSXX" and creates the loss matrix ror 
riscal year "ХХ". The resultaat natrix is "LXX" tor fiscal 
megr "хх". The function "MATRIX" could create a matrix or 
the fcllowing dimensicn 7x40x10x31x8 for 7 years, 40 MOS's, 
10 grades, 31 LoS's, and 8 loss types. However, due to 
limited workspace, the dimersion of 490x31x19 for ЧО MOS's, 
31 LCS'S, and 10 grades was commonly utiiized when the 
losses were aggregated. Each loss type could be examined 


Witn the sane dimensicred matrix. See Figure D. 15. 


V GETLOSS 
a THIS GETS THE L(XX) MATRIX FOR YEARS 


ΓΗ ΕΙΚ [OSS77 
( /8-МАТКІХ (05578 
МАК 5579 
(80-МАТКІХ (05530 
а(81--МАТКІХ (05581 
а| 82--МАТНБІХ 105532 
8(85-МАТКІХ 105583 


OO YD Q + YN — 
A A ر‎ A A en 


90 0 
| 
| 
| 
| 


Figure D. 14 APL Function GETLOSS. 


Een Groupe LOSS Fubctions 


There is a separate function to make a conbined 
Matrix tor tach of the MOS groups. Listed in Table XXIV are 
tae groups, the function used to create tne desired matrix, 


and the shape of the resultant matrix. See Figure D.16 ror 





the Ground Combat function. THe Otc: groups’ runctions 


were similar. 


| 
| 
| 


Ç GCLOSS 

[1] a THIS GETS THE LOSSES FOR GROUND COMBAT 
[2] e FOR 6 YEARS 77 - 82. 

[3] a SHAPE IS (MOS,YEAR,LOS,GRADE). 

[4] a CHANGE YEAR INDEX TO DESIRED NUMBER 
ی۷۷‎ ۶۶ + 70, 

[6] ٠ہ6ے]1:1::[“-77]4::[‎ 

15 СОИ] 
ЭСЕ] 

[9] GCL[1;4;;]-L80[ 4:;] 

[10] a ссі; 5|-081(4;;| 

[11] а ССЦ 16;|-182(4;;| 

БІРЛ ο ο L776; | 

ЕЕЕ) е672[65: 1] 

[i me c [2 s J=etL7ə9[6:] 

[15] 6CL[2:4::]-L80[6::] 


[16] a GCL[2;5;;] - LB1[ 6; 
[17] a GCL[2;6;;]-1L82[6; 


|] 
) 
ΠΟΤΕ] 
[19] GCL[3;:2;;]-L78[11;;] 
Оа зеби] 
[επ ο δα) 


[22] a GCL[3,5;;]-L8 1[11;:] 
[23] « 6CL[3:6;;]- L82(1 1;;] 


== A ne, EM Gee | — E MA ο, گلھت‎ πε... III σας... «ος. πας... πα. a θες... iii ¿at ӘШБ. дра, 
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Figure D.16 APL Function GCLOSS. 


These functions use the global variables "LX", 
ШО еасп "ХХ" Tiscal year, In order to get a Combat Service 
Support HOS loss matrix, the "LXX" matrix was used directly 
to remove tne desired comLination or celis. 
To create the loss Matrix used tor the 
Janes-Stein function the rollowing functions were used: 
КО ο ecrcates the "YC" matrix for Conmvat Support. 
ο στ στ οτε the "Y1G" matrix ror Ground Comtat. See 


Figure D.17 for an exaupie. 





e let K = nunber of feasible cells, i.e., sum of all Dij. 
e let Yij = number of leavers in cell (i,j). 

e Let t = 1,...,T; where T = number of years of data used 
to create the estimator. 

° let Nij = Inventory in ceil (153) Ξ Max 
((N(t) +N (t+1))/2,Y(t))- 

The following АР, functions were utilized in the 

James-Stein model. See Figure D.18 for the APL listings of 


the fcllowing functions. 
a. Function ARCSIN 


Let ISM = Nij ARCSIN PMMij. The function 
"ARCSIN" returns the inverse sine transformation for use 
when the success probability is estimated directly (e.g., by 
MINMAX). 


b. Function EINCONV 


Let Pij = Jij BINCONV Nij(t). The function 


"BINCCNV" inverts the inverse Sine transtormaticn. 
c. Function EINPREP 


Let ISij(t) = Yij(t) BINPREP Naij(t). The func- 
tion "BINPREP'" prepares the Freeman-Tukey version or the 


inverse sine transformation for binomial data. 
d. Function JAMES 


Ше II). DI] JAMES 21515 (t). The function 
"JAMES" returns a James-Stein estimator for the means of the 
means of the cells in the last two dimensions of the input 
Matrix while being screened by the incidence matrix of 


nonstructural zeroes. 
e. Function MINMAX 


Let PMMij = Yij(t) MINMAX Nij(t). Tne function 
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function was utilized in the validation procedures to 
compare the James-Stein estimator to the MINMAX estimator. 


Seer igure Da 19 for the AFL listing of this function. 


— — — TO, A CA GEE gs cece unmmq n x. —. — -. — — E — a Lnd 


V P<Y MINMAX N;YY 
[1] a RETURNS MINMAX EST'S FOR BINOMIAL 
[2] a INCLUDES REDUCTION OVER THE 1st RANK (TIME) OF THE 
[3] ^ INPUTS. NN (GLOBAL) IS CREATED FOR USE IN ARC SINE. 
[4] vYY--47Y 
[5] NN<—+/N 
[6] P-(0.5*(OzYY)xYY-NN*0.5) 2 12 NN«0O.5 


P COMER UNE — — — sh . — I MS aU ‘M I BS шщ сше шо ш a A <  ———— оаа Hs ERE — s 


Figure D.19 APL Function MINMAX. 
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4. Validation Functions 


= ce Se ss Qm sss ο... 


The following APL functions are an example of the 
functions utilized in the validation procedures specified in 
ο ντε 5. Each MOS group (Aviation, Coubat Support, Combat 
Service support, and Ground Combat) was unique in its set or 
fonctions. However, the metncdology was Similar 1n nan1pu- 


lating the data. 
a. Function BEFORE 


The function "BEFORE" performed the preliminary 
work for the parameters ror the James-Stein, MinLiaax, and 


Maximum Likelihood estimation processes. See Figure D.2). 
L. Function GETGCBI 


The functicn "GETGC81" derives the Normal values 
for the James-Stein, Minimax, and Maximum Likelihood estıma- 


COrs. See figure D.2l. 
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μονο A A ست ہے‎ шшш A A E es ge, A A a O rn E O nn A RD. ARS A ME | шыу, سوک ,ہے ھکوس‎ A AS «ңы KB ӘӘ аш... m 


— — r a r s s rÑ w — παι 


η; S 51 Y.NSNS8:YƏ.YS IOS. O8 
] a GETS THE NORMALS FOR ΤΗΕ JAMES-STEIN, MINMAX, 
| a AND MAXLIKE FOR THE GROUND COMBAT MOS FOR THE 
| a PROJECTED YEAR 1981. ALSO COMPUTES THE CHI 
] a SQUARE VALUES FOR VALIDATION PURPOSES. 
] 1-0 
| М5-М8165 
| М8-М8168 
ГІ 5202202: 
| ۲88 
0| 165-(Ү5 8۱۱۷۸۲۲ М5) 
1] 168-(Y8 BINPREP NS) 
2] NOR81G5+-1G5-JGS 
3] NOR8168-168-J68 
4] DD-DLT 
5) P-PG5 
[16] С5816-(2 6)р0 
ШОО ESE σης NS 
[18] 00-0LC 
[19] P-PGB8 
[20] CS81G[2; 1 2]-Y8 CHI N8 
[21] NMM8165-165-1M5 
[22] NMM8168-168-1M8 
[23] DD-DLT 
] P-PM5 
P l [i rr 3] y5 CHINS 
] 0D-DLC 
] P-PM8 
| С5816(2; 3 4]-Y8 CHI N8 
| МА8165-165-ІА5 
| МА8168-168-ІА8 
1 00-011 
] P-PA5 
ese Gili 5 6|-ү5 СНІ М5 
] DD«-DLC 
| Р-РАВ 
6581605 6J=Y8 CHI N8 
|! JAMES DF MINMAX DF MAXLIKE OF’ 
| 0-С5816 


Ç 
| 
| 
| 
| 
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Figure D.21 APL Function GETGCBI. 
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ТЕЗІ 15A.B C:DB:1S;J;P:MLE;AC;SSEM, SEM, SSEP:SEF;DIFP;DIFM.NUM 
Y [5 LEAVERS, SHAPE (2 YEARS, 31 LOS, MOS'S) 

D IS THE GLOBAL INCIDENT MATRIX FOR THE DESIRED GRADE 

N IS CENTRAL INVENTORY, SHAPE (2 YEARS, 31 LOS, MOS'S) 
IS-Y BINPREP N 

JeD JAMES IS 

P-(J BINCONV N)xO 

MLE- COT 1) 7 (NE Tz] (NC 1::]20))x0 

ocu an (( ii) i ((r/N)=0)))xD 

SSEM*- */, SEM - (AC -MLE)*2 


SSEP - */, SEP -(AC-P)*2 
DIFP-+/,(SEP<SEM) 
DIFM- +/,(SEM<SEP) 
NUMe+/,D 
IAS TESTO. ALL TRE CELLS IN THE FEASIBLE REGION’ 
' ARE ESTIMATED BY THE JAMES STEIN PROCESS' 
SEP = THE SQUARED DIFFERENCE OF THE JS PROJECTEO AND ACTUAL' 
' SEM = THE SQUARED DIFFERENCE OF THE MLE AND ACTUAL' 
' SSE(MLE) = THE SUM OF THE SQUARED DIFFERENCES OF SEM’ 
ΠΤ ΞΙΊ ΕΞ THE SQUARED DIFFERENCES OF SEP’ 
‘FEASIBLE CELLS SEP<SEM SEM<SEP SSE(MLE) SSE(P)' 
Α-' | 
Be' ' 
Ce 
0 ۲٦ 


B,(SNUM),A,(3DIFP),8,(+DIFM),C,(3SSEM),DD,($SSEP) 


SSEM e -/, SEM-(AC - (MLEx (AC 20)))*2 

SSEP e */,SEP - (AC -(Px(AC0)))*2 
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NUM+e-+/,(AC #0) 
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ARE REPRESENTED IN THE COMPARISON TEST °‏ 


πο ο τς SQUARED DIFFERENCE OF THE JS PROJECTED AND ACTUAL’ 


' SEM2 = THE SQUARED DIFFERENCE OF THE MLE AND ACTUAL’ 

' SSE2(MLE) = THE SUM OF THE SQUARED DIFFERENCES OF SEM2' 
' SSE2(P) = THE SUM OF THE SQUARED DIFFERENCES OF SEP2' 
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Figure D.24 APL Function TEST. 
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